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Abstract. We consider the free boundary problem for a layer of viscous, incompressible fluid 
in a uniform gravitational field, lying above a rigid bottom and below the atmosphere. For 
the "semi-small" initial data, we prove the zero surface tension limit of the problem within a 
local time interval. The unique local strong solution with surface tension is constructed as the 
limit of a sequence of approximate solutions to a special parabolic regularization. For the small 
initial data, we prove the global-in-time zero surface tension limit of the problem. 



^ ■ 1. Introduction 

^ I 1.1. Formulation of the equations in Eulerian coordinates. We will follow Guo and Tice 

[12, 13, 14] to formulate our problem. We consider a viscous, incompressible fluid evolving in a 



moving domain 

n{t) = {y e S X M I - b{yuy2) < ya < ^2, t)}. (1.1) 
Here we assume that either S = M^, or S = (LiT) x (L2T) for T = R/Z the usual 1— torus and 
Li,L2 > the periodicity lengths. The lower boundary of ^{t), denoted by S^, is assumed to 
be rigid and given, but the upper boundary, denoted by S(t), is a free surface that is the graph 



(1.2) 



in 

cm: 
< 

■ of the unknown function r? : S x — t- M. We assume that 

^— > ■ 

^' fo < 6_ < 6 G C°°(S) if S = (LiT) X (L2T 

[0 < 6_ < 6 G C°°(S) and 6 ^ 600 > as |(yi,y2)| ^00 if S 

for some constants 6_,6oo- The fluid is described by its velocity and pressure functions, which 
'i> [ are given for each t > by u{-,t) : Q,{t) — )• and p{-,t) : Q{t) — )• M, respectively. For each 

^ ■ t > we require that {u,p,r]) satisfy the following equations 

dtu + u ■ Vu + Vp — fj,Au = in Q,{t) 

divM = in ^}(t) 

{pi - fiB{u))u = gr]v - aHv on T.{t) (1.3) 

dtV = U3 - uidy^r] - U2dy2r] on $](*) 

M = on Eft 



for / the 3x3 identity matrix, (Bti)jj = diUj + djUi the symmetric gradient of u, /i > the 
^ ■ viscosity, and g > the strength of gravity. The tensor pi — jjl]){u) is known as the viscous 

5^ . stress tensor. We take cr > to be the constant coefficient of surface tension, and it should 

be understood that when a = it means the case without surface tension. In this paper, we 
let V,, denote the horizontal gradient, div^, denote the horizontal divergence and A,, denote the 
horizontal Laplace operator. Then the outward-pointing unit normal on S(t), v, is given by 

(-V*r?,l) 



V1 + |V*?7P' 

and H, twice the mean curvature of the surface S(t), is given by the formula 



(1.4) 



iJ = div, [ ^ ) . (1.5) 
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The fourth equation in (1.3) is called the kinematic boundary condition which implies that the 
free surface is advected with the fluid. Note that in (1.3) we have shifted the gravitational 
forcing to the boundary and eliminated the constant atmospheric pressure, Patm, in the usual 
way by adjusting the actual pressure p according to p = p + gya — Patm- Without loss of 
generality, we may assume that ji = g = 1. For a more physical description of the equations 
(1.3) and the boundary conditions in (1.3), we refer to [26]. 

To complete the formulation of the problem, we must specify the initial conditions. We 
suppose that the initial surface S(0) is given by the graph of the function 'q{Q) = rjo : ^ M., 
which yields the initial domain 17(0) on which we specify the initial data for the velocity, 
m(0) = Mo '■ 17(0) — )■ M^. We will assume that % > — 6 on S and that (tiO)%) satisfy certain 
compatibility conditions, which we will describe in detail later. When S = (LiT) x (L2T) we 
shall refer to the problem as the "periodic case", and when E = we shall refer to it as the 
"non-periodic case" . In this paper, for the local-in-time analysis we will consider both the non- 
periodic and periodic cases simultaneously. When different analysis is needed for each case, we 
will indicate so. Otherwise, the argument we write works in both cases. For the global-in-time 
analysis we will only consider the periodic case for simplify. 

1.2. Geometric form of the equations. The movement of the free boundary E(t) and the 
subsequent change of the domain ^{t) create numerous mathematical difficulties. To circumvent 
these, as usual, we will transform the free boundary problem under consideration to a problem 
with a fixed domain and fixed boundary. We will not use a Lagrangian coordinate transformation 
as in [5, 20], but rather a fiattening transformation introduced by Beale in [6]. To this end, we 
consider the fixed equilibrium domain 

Q := {x G S X M I - 6(xi, X2) < X3 < 0} (1.6) 

for which we will write the coordinates as x G 17. We will think of S as the upper boundary of 
O, and we continue to view 77 as a function on S x R+. We then define 

fj := Vrj = harmonic extension of rj into the lower half space, (1.7) 

where Vrj is defined by (2.1) when S = and by (2.2) when E = (LiT) x (L2T). The harmonic 
extension fj allows us to fiatten the coordinate domain via the mapping 

Qb {xi,X2,X3 + 7?(x, t){l + X3/b{xi,X2))) := $(x, t) = (yi, 2/2, ys) G ^{t). (1.8) 

Note that $(E,t) = E(t) and $(-,t)|s(, = Id-£^, i.e. $ maps S to the free surface and keeps the 
lower surface fixed. We have 

/I 0\ /I -AK\ 

V$ = 1 and ^ := (V$~^)^ =01 -BK (1.9) 

\A B j) \0 K J 

for 

A = difjb — {x3f]dib)/b'^, B = d2fjb — {x3f}d2b)/b^ , 

J =l + fj/b + d3fjb, K = J-\ (1.10) 
6= (I + X3/6). 

Here J = det V$ is the Jacobian of the coordinate transformation. 

If rj is sufficiently small (in an appropriate Sobolev space), then the mapping $ is a diffeo- 
morphism. This allows us to transform the problem to one on the fixed spatial domain Q for 
t > 0. In the new coordinates, the system (1.3) becomes 



dfU — dtrjbKd^u + u ■ V^u — A^n -I- = 


in 


div^ u = 


in Q, 


Sa{p,u)M = r]M -aHM 


on E 


dtrj = u ■ M 


on E 


u = 


on E;, 


Xu,ri) |t=o= {uo,m)- 





(1.11) 
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Here we have written the differential operators V^, div_4, and with their actions given by 
(V^/)i := Aijdjf, div^X := AijdjXi, and A^/ := div^V^/ for appropriate / and X. We 
have also written M := (—diij, —927/, 1) for the non-unit normal to and we have written 

Sa{p-i '■= pI—^AU for the stress tensor, where (]D)^u)ij := AikdkUj+AjkdkUi is the symmetric 
^—gradient. Note that if we extend div_4 to act on symmetric tensors in the natural way, then 
div^ S_A.ip, u) = V_AP — ^A'^ vector fields satisfying div^ u = 0. 

Recall that A is determined by rj through the relation (1.9). This means that all of the differ- 
ential operators in (1-11) are connected to r], and hence to the geometry of the free surface. This 
geometric structure is essential to our analysis, as it allows us to control high-order derivatives 
that would otherwise be out of reach. This was well explained in [12, 13, 14]. 

1.3. Main results. Free boundary problems in fluid mechanics have been studied by many 
authors in many different contexts. Here we will mention only the work most relevant to our 
present setting. For a more thorough review of the literature, we refer to the review paper of 
Shibata and Shimizu [19] and the references therein. 

The problem of describing the motion of an isolated mass of viscous fluid bounded by a free 
boundary was studied by Solonnikov in a series of papers; for instance, [20] considers the problem 
without surface tension, while [21, 22] concern the problem with surface tension. Solonnikov's 
technique for proving the well-posedness in these papers did not rely on energy methods, but 
rather on Holder estimates in the case without surface tension and Fourier-Laplace transform 
methods in the case with surface tension. A local well-posedness for the problem with surface 
tension that based on the energy method was developed by Coutand and Shkoller [9] . 

The viscous surface wave problem in our present setting that describes the motion of a layer 
of viscous fluid lying above a fixed bottom has attracted the attention of many mathematicians 
since the pioneering work of Beale [5]. For the problem without surface tension, Beale [5] 
proved the local well-posedness in the Sobolev spaces and Sylvester [23] studied the global well- 
posedness by using Beale's method. For the periodic case with a completely flat fixed bottom, 
Hataya [16] obtained the global existence of small solutions with an algebraic decay rate in time. 
Recently, Guo and Tice [12, 13, 14] developed a new two-tier energy method to prove the global 
well-posedness and decay of this problem; they proved that the solution for the non-periodic case 
decays to the equilibrium at an algebraic rate, while the solution for the periodic case decays at 
an almost exponential rate. For the problem with surface tension, Beale [6] proved the global 
well-posedness in the Sobolev spaces, while Allain [1] obtained a local existence theorem in two 
dimension using a different method. Bac [4] showed the global solvability in Sobolev spaces 
via the energy method. Beale and Nishida [7] showed that the solution for the non-periodic 
case obtained in [6] decays at an optimal algebraic rate. For the periodic case with a flat fixed 
bottom, Nishida, Teramoto and Yoshihara [18] showed the global existence of small solutions 
with an exponential decay rate. Tani [24] and Tani and Tanaka [25] also discussed the solvability 
of the problem with or without surface tension within Beale-Solonnikov's functional framework. 

Formally, when the surface tension coefficients a tend to zero the solutions of the problem 
(1.11) with (T > should converge to the limit that solves the problem without surface tension. 
However, there is no general theory that guarantees such asymptotic limit; it depends on whether 
one can develop an appropriate well-posedness theory for the problem which is unified in the 
surface tension. We note that the well-posedness theories for the problem with surface tension 
mentioned above highly depend on the surface tension, so it can not be applied to show the zero 
surface tension limit. Motivated by the recent works of Guo and Tice [12, 13, 14], in this paper 
we shall develop a unified framework of nonlinear energy method of proving the well-posedness 
of (1.11) with and without surface tension. Hence, it will allow us to rigorously justify the zero 
surface tension limit in (1.11). Moreover, for the small initial data this limit is showed to be 
global in time. We remark that the problem of the zero surface tension limit was considered in 
the context of the water wave problem by Ambrose and Masmoudi [2, 3] and the Stefan problem 
by Hadzic and Shkoller [15], and both of them are restricted within a local time interval. 

Before we state our results, let us now mention the issue of compatibility conditions for the 
initial data (uq, rjo). We will work in a high-regularity context, essentially with regularity up to 
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2N temporal derivatives. This requires us to use uq and rjo to construct the initial data dfu{0) 
and 9^77(0) for j = 1, . . . , 2N and dlp{0) for j = 0, . . . , 2N — 1. These other data must then 
satisfy various conditions (essentially what one gets by applying df to (1.11) and then setting 
t = 0), which in turn require uq and 770 to satisfy 2N compatibility conditions. We describe 
these conditions in detail in Section 3.3.1 and state them explicitly in (3.83). 

To state out our results, we now define some quantities. We shall recall from Section 1.4 our 
notation for Sobolev spaces and norms. For a generic integer n > 3, we define the energy as 

n „ n— 1 n n+1 ^ 

(1.12) 



i=o 



j=0 



2n-2j-l 



2n-2j+3/2 



and the corresponding dissipation as 

2 



j=0 



2n-2j+l 



+ 



n-1 
j=0 



2n-2j 
n+1 



l^ll2n+3/2 + f^^ll^«^ll 2n+l/2 



(1.13) 



+ Mln-l/2 + \\dtV\\ln-l/2 + T.\n^ 



j=2 



2n-2j+5/2 



For our use, we will take both n = 2N and n = N + 2. We define 

£S ■■= \K\\In + ll%"ll4iV+i/2 + Iko ll4Ar+i ■ (1-14) 

Note that we allow for a = in the definitions of energy functionals throughout the paper. We 
also define 

-^27V := ||??|l4iV+l/2- (1-15) 

Wc first state the local-in-timc results for the "semi-small" initial data; we take either S = 
or S = (LiT) X (L2T). For any integer > 3, we define 



JC^N(t)-= sup ^2^(0+ / ^^Nir)dr + 

0<r<t Jo 



,^,,+ sup ^2Jv(r), 

{XtY 0<r<t 



(1.16) 



where we have referred to the space Xt defined later in (2.21). 



Theorem 1.1. Let N > 3 be an integer. Assume that the initial data Uq and tjq satisfy the 
bound £q < 00 as well as the (2iV)*'^ compatibility conditions (3.83) with Q < a < I. There 
exists a universal eo > such that if ||??o ll4Af-i/2 — ^o/2, then there exists a Tq = Tq{£q) > 
so that there exists a unique solution [u" ^p" ,r]") to (1.11) with initial data {uQ,r]Q) on the time 
interval [0, Tq]. The solution [u" ,p" ^rj^) obeys the estimates 



/C^iv(«",K,»?")(ro) < P{£^) and sup ||r?'^(t)||:^_i/2 < £0 

0<t<To ' 



(1.17) 



for a universal polynomial P with Pifi) = 0. 

As a consequence, if as a ^ 0, Uq — t- uq in H^'^ (ft), tjq 
in H^^~^^{T,), then {u" ^p'' ,rj") converges to the limit {u,p,r]) which is the unique solution to 
(1.11) for (7 = with initial data {uo,rjo) on [0,To]. The limit {u,p,r]) obeys the estimates 



r/o in i74^+i/2(S) and ^ril 



K.\^(u,p,T]){j:o) < P{£o) and sup 

0<*<To 



4Af-l/2 



< £0 



(1.18) 



We now state the global- in-time results for the small initial data; we take S = {LiT) x (L2T). 
We assume further that the initial surface function satisfies the "zero-average" condition 



^1-^2 Jy. 

For any integer iV > 3, we define 

g^N{t):= sup £SN{r)+ fv^N{r)dr+ sup {1 + r)^^-^£^^^{r) + sup ^ 

0<r-<t Jo 0<r<t 0<r<t U 



J^2N{r) 



+ r) 



(1.19) 



(1.20) 
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Theorem 1.2. Let N > 3 be an integer. Assum,e that the initial data Uq and tjq satisfy the 
{2NY^ compatibility conditions (3.83) with < a < 1, and that 7]^ satisfies the zero average 
condition (1.19). There exists a universal 6o > such that if £q < ^O; then there exists a unique 
solution (vF .,jf ,rf^ to (1.11) with initial data (ugj^o) ^'^ [0, oo). The solution {vf^ ,r]'^) obeys 
the estimate 

afivK,P",^")(oo)<Cfo" 

for a universal constant C > 0. 

As a consequence, if as a ^ 0, Uq — t- uq in H'^^ {Vl), rj^ — t- r/o in i7'^^+^/^(S) and y/^rjQ — 
in H^^~^^{Ti), then {u" ^p'' ,r]") converges to the limit {u,p,ri) which is the unique solution to 
(1.11) for (7 = with initial data (no,?7o) on [0, oo). The limit {u,p,r]) obeys the estimate 

Q%{u,p,r]){^)<C£o. (1.22) 

Remark 1.3. In Theorem 1.1 (resp. Theorem 1.2) the convergence of [u'^ ,p" ^rj") to {u,p,r]) 
is weak in the norms in the definition of lC2jq (resp. G2n) '^^^ strong in the norms of any 
functional space that /C^jy (resp. G2n) compactly embed into. The convergence is more that 
sufficient for us to pass to the limit as a in {1.11). 

Remark 1.4. Note that in Theorem 1.1 we only assume that the i74iv-i/2 

770 is small, 

but the higher-order norm of rjQ and the norm of uq can be arbitrarily large. In this sense we 
may regard the initial data as to be "semi-small". Theorems 1.1-1.2 contain the local well- 
posedness of the viscous surface wave problem (1.11) both with and without surface tension for 
the "semi-small" initial data and global well-posedness for the small initial data. Moreover, 
the solution without surface tension decays at an almost exponential rate and the solution with 
surface tension decays at an exponential rate (cf. Theorem 5.14)- Since rj is such that the 
mapping ^(-jt), defined by (1.8), is a diffeomorphism (C^^~^/^ for a > and (74JV-2 j^^, 
a = 0) for each t > 0, we may change coordinates to y E fl{t) to produce solutions to (1.3). 

We now explain the strategy of proving Theorems 1.1 and 1.2. After collecting some prelimi- 
naries in Section 2, we start to prove Theorem 1.1 of the local-in-time results. The key step is to 
establish a unified local well-posedness theory for the problem (1.11) with and without surface 
tension. We manage to first produce the local unique strong solutions to (1.11) for each fixed 
a > 0. In the framework of the energy method, it is more natural to use the geometric formula- 
tion (1.11) rather than the perturbed form as well explained in Guo and Ticc [12, 13, 14]. For 
the local well-posedness without surface tension, Guo and Tice [12] used an iteration scheme 
based on the solvability of the two linear problems which are separated from (1.11), that is, the 
linear ^-Navier-Stokes equations for {u, p) , where A and M are given (in terms of 77) , 

' dtu — A_4ti + V^p = in Q 

div^ u = in J7 

' Sa{p,u)M = F^ onS 

u = on Eft, 

and the linear transport equation for 77, where u is given, 

dtri + u- V^T] = U3 on S. (1.24) 

Hereafter u ■ Vi^ri := uidirj + U2d2ri, and etc. A subtle point for the iteration in [12] could be 
useful is that for the nonlinear problem (1.11) without surface tension, F"^ = rjM (~ Vf] ~ Vu); 
this would imply that the regularity of rj obtaining from (1.24) {rj ~ /g ^^3) can control the 
nonlinear term F"^. However, for the problem (1.11) with surface tension, F^ = r]J\f — aHJ\f 
(~ V^?7 ^ Vu); then the regularity of rj obtaining from (1.24) is insufficient to control the 
nonlinear term F^ . Hence, we can not apply the iteration scheme of [12] directly. One may 
then instead use an iteration scheme based on the solvability of the following linear problem for 



1.21 



(1.23) 
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{u,p,r]), where A and Af are given (in terms of, say, Q, 




(1.25) 



We may remark here that wc could solve the problem (1.25) in a weak sense by using a time- 
dependent Galerkin method similarly as [12] (restricted to the div^-free vectors), however, the 
standard way of improving the regularity by using the horizontal difference quotient as in [9, 4] 
would fail to improve the regularity of the weak solution due to that the horizontal difference 
quotient does not commute with the projection along the normal J\f. 

Our way of getting around the difficulty explained above is to introduce an artificial viscosity 
term in the transport equation, that is, we consider the following regularized K-problem, where 
K > is the artificial viscosity coefficient. 



Note that we have also introduced along the so-called compensator function ^ in the transport 
equation owing to the issue of compatibility conditions for the initial data. The function ^ 
is determined by the initial data, and we postpone its construction in Section 3.3.2. By the 
introduction of such then at time t = 0, we essentially add nothing on the boundary, and 
hence the compatibility conditions to (1.26) are same as (1.11). This allows us that the initial 
data {uo,r]o) for the original problem (1.11) can be exactly taken as the one for the regularized 
problem (1.26); otherwise, even though (ncryo) satisfy the compatibility conditions for (1.11), 
in general they do not satisfy the corresponding compatibility conditions for (1.26). One may 
refer to [8] for more discussions about the issue of compatibility conditions. It is worth noting 
that we employ a different construction of the compensator * from [8]. The way of [8] is to 
define the Taylor polynomial of — A*77(t) at t = to be the compensator which requires more 
regularity on the initial data; hence a very complicated argument of regularizing the initial data 
that ensures the compatibility conditions needs to be carried out. Our construction avoids this. 

By adding this artificial viscosity term, the regularity of rj obtaining from the fourth equation 
of (1.26) (A*r/ ~ U3) can control the nonlinear term rjj\f —ctHM . This will allow us to modify the 
iteration scheme of [12] to establish a local well-posedness of the ^-problem (1.26) in Section 3 
based on the solvability of the linear ^-Navier-Stokes equations (1.23) and the following linear 
regularized surface problem 



Having obtained the unique strong solutions {u'^,p'^, rj'^) to (1.26) on a local time interval [0, Tq] 
for each k > (cf. Theorem 3.13), in Section 4 we then derive the K-independent estimates for 
{u'^,p'^,r)'^) which allow us to pass to the limit as k ^ in (1.26) to produce the unique strong 
solution ('u^,p'^,77^) to (1.11) on a local time interval [0, Tq'] for each fixed a > (cf. Theorem 
4.9). We next continue to derive the fi-independent estimates for [u" ^p" ,r]") which allow us to 
pass to the limit as a — )• in (1.11) within a local time interval [0, Tq]. In order to derive the 
(T-independent estimates (and the derivation of the K-independent estimates is similar), it is a 
very key point to design the special quantity (defined by (4.96)) which allows us to establish 
the following estimate 



< 




(1.26) 



dtT] - kA^ij = + on E. 



(1.27) 




(1.28) 
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Since (^0''7o) ~^ ^o('"0)%) as cr — )• 0, a standard continuity argument basing on (1.28) then 
implies that the solutions {u'^ ,p'^ ,rj"') to (1.11) all exist on a u-independent time interval [0, Tq] 
and satisfy < 2P{£q ). We can then improve this estimate to be (1.17) by improving the 
estimates of 77. Note that we will refine many estimates of [12] so that Theorem 1.1 can be 
established for our "semi-small" initial data which is not necessary to be genuinely small. This 
completes the proof of Theorem 1.1. 

We then proceed to prove Theorem 1.2 of the global-in-time results in Section 5, and the key 
step is to derive the global-in-time cr-independent estimates for the local solutions implicitly 
contained in Theorem 1.1. We will employ the two-tier energy method recently developed by 
Guo and Tice [13, 14]. In this paper, for simplify we choose to consider the "periodic case" 
but we could treat the "non-periodic case" by using the arguments of [14]. For the sake of 
completeness, here we will briefly sketch the two-tier energy method. Basing on the natural 
energy evolution and making full use of the structure of the problem (1.11), we may prove the 
energy estimate, where /C < £"^+2^ 

£^N{t) + f T^2N ^ ^2^(0) + /CJ-2iv + /* V^2iv^-^2Ar + ICT2N ■ (1-29) 

JQ Jo 

The derivation of (1.29) is divided into two steps. First to use the energy evolution, we can only 
consider the "horizontal" energy and dissipation with localization, say, £^ and P„ (cf. Section 
5.2 for the definitions), and we may prove 

4V(i) + /* ^2iV < f 2^(0) + (f 2V(i))'/' + /* V^'^2N + V^27V^-^2iV + sV^^. (1.30) 

Jo Jo 
Second we use the structure of (1.11) to show the comparison estimates that 

< ^2N + ^^2N. and V^^ < V2N + ICT2N- (1.31) 

It is easy to show that fgAr is comparable to stated in (1.31) by applying the Stokes 

elliptic theory of Lemma 2.15. However, due to the lack of t] terms in T>2n, we can not compare 
the dissipations by using this technique. To overcome this difficulty, Guo and Tice [13, 14] 
instead use the structure of the equations as well as the equation for the vorticity (from which 
the pressure and rj can effectively be eliminated) to derive various estimates for u and Vp 
without reference to 77. Then a bootstrapping procedure is employed to obtain estimates for 
rj and p (not just its gradient). In this paper, we use an alternative procedure which seems 
much more simpler. More precisely, since T>2N controls horizontal derivatives, we can use it to 
gain the regularity of w on S. The idea is then to apply the Stokes elliptic regularity theory of 
Lemma 2.16 to deduce the desired estimates of u. Then we get the desired estimates of p and 
r] by using the equations directly. Then (1.29) follows by (1.30) and (1.31). 

In order to close the estimates (1.29), we appeal to the control of J-2N- Note that for a > 0, 
we can bound by J^2N < '^ '^'^2N- Hence, for each fixed a > 0, from (1.29) we may prove the 
global well-posedness and exponential decay of (1.11) for the initial data sufHciently small with 
respect to a (cf. Theorem 5.14). However, to prove the zero surface tension limit, we need to 
derive the a-independent estimates for the solutions. The only way to estimate F2N then is 
through the kinematic transport equation for 77, and we may derive 

J'2N{t) < J^2n{0) + t [ V^2N- (1-32) 

Jo 

This estimate allows F2N to grow linearly in time. So to balance such growth of J^27Vj we must 
derive a strong decay of /C < £n+2- This is accomplished by the following differential inequality 

j/n^2 + '^n+2<^- (1-33) 

Indeed, we interpolate to have f^r+a ^ {Vli^2)^^^-^^l^'^^-'^\£^j^f/^*^-'^\ so we can deduce 
from (1.33) that 

(1 + t)4^-«fXr-H2(i) < ^2^(0) +^2Ar(0). (1.34) 
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Then (1.29), (1.32) and (1.34) imply the desired bound Q^j^{t) < afjv(O). This a-independent 
estimates will then allow us to pass to the limit as cr ^ in (1.11) on the whole interval [0, oo). 
Hence we can complete the proof of Theorem 1.2. 

1.4. Notation convention. Wc shall also employ the notation convention from [12, 13, 14]. 

We write i?'=(Jl) with A; > and H'{T,) with s G M for the usual Sobolev spaces. We also 
typically write = L^. To avoid notational clutter, we will avoid writing H''{Q,) or ii"^(E) 
in our norms and typically write only Since wc will do this for functions defined on both 
0, and S, this presents some ambiguity. Wc avoid this by adopting two conventions. First, we 
assume that functions have natural spaces on which they "live." For example, the functions u, 
p, and fj live on J7, while 77 itself lives on S. As we proceed in our analysis, we will introduce 
various auxiliary functions; the spaces they live on will always be clear from the context. Second, 
whenever the norm of a function is computed on a space different from the one in which it lives, 
we will explicitly write the space. This typically arises when computing norms of traces onto S 
of functions that live on Jl. We use denote the norm of the space 1^(0, T;X). 

Wc write N = {0, 1,2,...} for the collection of non-negative integers. When using space-time 
differential multi-indices, we will write N^"*""^ = {a = {ao,ai, . . . ,am)} to emphasize that the 
0— index term is related to temporal derivatives. For just spatial derivatives we write N"*. For 
a € N^"*"™ we write 9" = 9"° 9"^ • • • d^. We define the parabolic counting of such multi-indices 
by writing \a\ = 2a;o + ai + • • • + a„i. We will write V*/ for the horizontal gradient of /, i.e. 
V*/ = difei + 82/62, while V/ will denote the usual full gradient. 

For a given norm ||-|| and integers A; > m > 0, we introduce the following notation for sums 
of spatial derivatives: 



E 

in<\a\<k 



5"/f and 



^ mJ 



(1.35) 



m<\a\<k 



For space-time derivatives we add bars to our notation: 



V ^ f 



aeNi+2 
m<|a|<fc 



^ mJ 



(1.36) 



m<.\a\<k 



We allow for composition of derivatives in this counting scheme in a natural way; for example. 



*m+'l 



f 



(1.37) 



Throughout the paper we assume that iV > 3 is an integer. We will employ the Einstein 
convention of summing over repeated indices. Throughout the paper C > will denote a generic 
constant that does not depend on the data, the surface tension coefficient a and the artificial 
viscosity coefficient k, but can depend on the other parameters of the problem, N > 3 and Q. 
We refer to such constants as "universal". Wc also use P to denote for a universal positive 
polynomial with P(0) = 0. Such constants and polynomials are allowed to change from line to 
line. When they depend on a quantity z we will write Cz and Pz to indicate this, or we will 
explicitly point out along the context. We will employ the notation a < 6 to mean that a < Cb 
for a universal constant C > 0. 



2. Preliminaries 

2.1. Analytic tools. 

2.1.1. Poisson integral. For S = R^, we define the Poisson integral in 



X (-00,0) by 



7M2 



(2.1) 
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where we have written /(^) the Fourier transform of / in M^. For S = {LiT) x (L2T), we define 
the Poisson integral in (LiT) x (L2T) x (— oo,0) by 

rf{x)= ^2^in-x'^2^\n\x3j^y^^ (2.2) 

where for n G {L^^'L) x (L^^Z) we have written 

fin) = / f{x'f-——dx'. (2.3) 

It is weh known that V : iJ*(E) i?*+i/2(0_) for 0_ = E x (-00, 0) is a bounded hnear 
operator for s > 0. However, if restricted to the domain we have the following improvements. 

Lemma 2.1. Let Vf he the Poisson integral of a function f that is either in ij*(S) or ii"5~^/^(S) 
for q eN, where is the usual homogeneous Sobolev space of order s. Then 

W^'-PfWl < and WV^VfWl < ||/||^,(^) . (2.4) 

Proof For the proof, we refer to Lemma A. 7 of [12] for the case S = and Lemma A. 9 of [12] 
for the case S = (LiT) x (L2T). □ 

2.1.2. Some inequalities. We need some estimates of the product of functions in Sobolev spaces. 

Lemma 2.2. Let the domain of function spaces be either S or Q,. Let < r < si < S2 be such 
that S2 > r + n/2, then the following hold when the right hand side is finite: 

WfgWHr < WfWHn Mh^. (2.5) 

and 

WfgWn-n < \\f\\H-r \\9\\h^2 ■ (2.6) 
Proof. Sec Lemma A.l of [12]. □ 
Remark 2.3. We will use a slight variant from (2.6) when the domain is for s > 5/2, 

ll/5ll(oHi(Q))* ~ \\f\\{om(n))* llfl'llH^(n)' (2.7) 
where the space oH^{Q,) is defined in (2.15). 

We then record some Poincare-type inequalities for the domain 
Lemma 2.4. Let f € H^{^), then 

ll/llL2(n)^ll/IL2(s) + l|93/IL2(n). (2.8) 
Let f G H^(Q,) be so that f = on T^b, then 

II/IIl2(Q) ~ 

(2.9) 

Proof See Lemmas A. 12 and A. 14 of [12]. □ 

We will need the following version of Korn's inequality. 
Lemma 2.5. It holds that ||n||^ < ||Bu||q for all u G iJ^(0;R^) so that u = on S^. 
Proof. See Lemma 2.3 of [6]. □ 

2.1.3. Continuity and temporal derivatives. In our local wcll-poscdncss arguments it is crucial to 
use the following interpolation result, which affords us control of the L°°H^ norm of a function 
/, given that wc control / in and dtf in L'^H^~'^. 

Lemma 2.6. Let the domain of function spaces he either S or Q. Suppose that f G -^^^([0, T]; i?*^ ) 
and dtf e L'^{[0,T]-H'^) for si> S2. Let s = {si + S2)/2. Then f e C^{[0,T]- H') (after pos- 
sibly being redefined on a set of measure ), and 

ii/iiioc^. < 11/(0)11^. + mw. + mwW^ ■ (2.10) 

Proof. See the proof of Lemma A. 4 of [12]. □ 
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2.1.4. Extension results. In our Local well-posedness arguments we need to be able to take the 
initial data (^u(O), j = 0,. . . , 2N and extend it to a function u° satisfying certain estimates. 

Lemma 2.7. Suppose that c^'u(O) G //^^^^■'(O) for j = 0,...,2N. Then there exists an 
extension u^, achieving the initial data, so that 

d{u^ e L\[o,ooy,H^^-^^+\n))nL^{[o,ooy,H'^^-^^{n)) 



for J = 0, . . . , 2N. Moreover, 

2N 



j=0 



u 



+ 



LooH4N-2j 



2N 



j=0 



2 

H4N-2j 



L2H4.N-2j+l 

Proof. See the proof of Lemma A. 5 of [12]. 

A similar extension result holds for c^?7(0), j = 0, . . . , 2N. 

Lemma 2.8. Suppose that ^?7(0) G H'^^-'^^+^{E) for j = 0,...,2N. Then there exists 
extension rf , achieving the initial data, so that 



(2.11) 

(2.12) 
□ 



an 



for j = 0, . . . , 2N. Moreover, 

2N 



l2jj4N-2j+2 



+ 



J^x, ]-[4N -2j + 



2N 

j=0 



fj4N-2j+l 



(2.13) 

(2.14) 

□ 



Proof. The proof is same as Lemma 2.7. 
2.2. Time-dependent functional setting. 

2.2.1. Time- dependent function spaces. We begin with introducing some function spaces. Define 

{ue H\n)\u\s,=0}, 

{u G H^{n) I u\s = 0}, and (2.15) 
{u G oH^in) I div'u = 0}, 



with the obvious restriction that the last space is for vector-valued functions only. 

For our time-dependent function spaces we will consider rj as given with A, J, etc determined 
by r) via (1.10). We define a time-dependent inner-product on = by introducing 



iu,v)y^O ■= (U- V)j{t) 



with corresponding norm ||?x||-^o := -s/ {u, u)-^o. Then we write n^{t) := 
larly, we define a time-dependent inner-product on qH^{Q,) according to 



{u,v)-^i := (B^(t)U : B^(t)u) J(t), 



(2.16) 

l^ll^o < oo}. Simi- 
(2.17) 



and we define the corresponding norm by = \/ {u, u)yr . Then we define 

'M}[t) := {u I ||n||^i < oo,u|s(, = 0} and X{t) := {u G ^^(t) | div_4(t) u = 0}. (2.18) 

We will also need the orthogonal decomposition 'H^{t) = y{t) © y{t)^, where 

y{t)^ ■.= {VA{t)V\v(^''H\^)}. (2.19) 

In our use of these norms and spaces, we will often drop the it) when there is no potential for 
confusion. Finally, for T > and = 0, 1, we define inner- products on L'^{[Q,T]-H^{Vl)) by 



{u 



^'")h% ■= I {u{t),v{t))^kdt. 
^ Jo 



(2.20) 
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Write II "ull^fc for the corresponding norms and Hj, for the corresponding spaces. We define the 
subspace of div^i-free vector fields as 

At := {uenrl divA{t) u{t) = for a.e. t e [0,T]}. (2.21) 

Under a smallness assumption on r/, we can have that the spaces H'^(t) and Hj' have the 
same topology as H'' and LF'H^ , respectively. 

Lemma 2.9. There exists a universal eQ> Q so that «/supo<«T ||^(*)||3 < ^O) then 

A= Mk < MIh'' ^ ^ ll^llfe (2-22) 
v2 

for k = 0,1 and for all t G [0, T]. As a consequence, for k = 0,1, 

\\u\\L2Hk < ||u||-^fc < V2\\u\\L2Hk . (2.23) 

Proof See Lemma 2.1 of [12]. □ 

Remark 2.10. Throughout the rest of this paper, we will assume the assumption of Lemma 2.9 
is satisfied so that (2.22) -(2.23) hold. 

The following lemma provides the H~^/^ boundary estimates of u knowing u, div^ u G 7i^{t). 
Lemma 2.11. Ifue VP{t) and divert G T-L^{t), then 

ll^--^ll//-i/2(s) + 11^- ^^11^-1/2(26) ~ + l|div^u||^o . (2.24) 

Proof. See Lemma 3.3 of of [12]. □ 

The following lemma records the differentiability of norms in the time-dependent spaces. 

Lemma 2.12. Suppose that u G dfU G {H^)* ■ Then the mapping t i->- ||u(t)||^o(t) is 
absolutely continuous, and 

j^\\u{t)f^o =2{dtu{t),u{t))^^^y + J^\u{t)f dtJ{t) (2.25) 

for a.e. t G [0,T]. Moreover, u G C^i[0,T]; H^{n)). If v G V.^,, dtv G {n\,)* as well, then 

- {u{t),v{t))^o = {dtu{t),v{t))(^^^^. + {dtvit),u{t))^^r^. + j^u{t) ■ v{t)dtJ{t). (2.26) 

A similar result holds for u G Xt with dfU G (Xt)*- 

Proof See Lemma 2.4 of [12]. □ 

It is more important to show that the space oH^{Q,) is related to the space X{t). To this 
end, we define the matrix 

/ K 0\ 

M := M{t) = KV^ = K • (2-27) 

\AK BK l) 

Note that M is invertible, and = JA^ ■ The matrix M{t) induces a linear operator 

Mt ■ u M.t{u) = M{t)u, the properties of which are recorded in the following. 

Proposition 2.13. For each t G [0, T], Ait is a hounded, linear isomorphism: from H^{Q) to 
H''{n) for k = 0,l, 2; from L'^{n) to VP{t); from oH^{n) to V}{t); and from oHl{n) to X{t). 
In each case the norms of the operators A4t,A4'i~^ are bounded by a constant times 1 + ||'7(i)||9/2- 
Moreover, the mapping Ai given by A4u(t) := A4tu{t) is a bounded, linear isomorphism: 
from L^{[0,T];H''{n)) to L'^{[0,T]; H''{n)) for k = 0,1,2; from L^{[0,T]; H^\n)) to n^; from 
L'^{[0,T];oH^{n)) to H^; and from L'^{[0,T];oH^{^)) to Xt- In each case, the norms of the 
operators M. and Ad"^ are hounded hy a constant times the sum 1 + supq<«t ||'7(i) 119/2- 
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Proof. We note that to show that M.t map div-free vector fields to div^-free vector fields, the 
following will be used: since J / 0, 

p = diYA v^Jp = JdivA V = dj{JAijVi) = dj{JA^v)j = dj{M-^v)j = div(M-^u). (2.28) 

See Proposition 2.5 of [12] for the proof of this proposition. □ 

2.2.2. Pressure as a Lagrange multiplier. The following proposition allows us to introduce the 
pressure as a Lagrange multiplier. 

Proposition 2.14. If At G {'H^{t))* is such that At{v) = for all v G X{t), then there exists 
a unique p{t) G IHPit) so that 

{p{t),dWAv)-^o = At{v) for all v G H^{t) (2.29) 

and \\p{t)\\j^o < (1 + lk(t)||9/2) ||At||(^i(t))*. 

If A E (H^)* is such that A[v) = for all v G Xt, then there exists a unique p G so that 

{p, divA = A{v) for all v eHt (2.30) 

and IIpII^o < (1 + supo<<<T h(t)||9/2) ||A||(^i^). . 

Proof See Proposition 2.9 of [12]. □ 
2.3. Elliptic estimates. 

2.3.1. Two classical Stokes problems. We recall the elliptic estimates for two classical Stokes 
problems with different boundary conditions. 

Lemma 2.15. Let r > 2. Suppose that f G H'-^in),he H'-^in) and G //''-3/2(S). Then 
there exists unique u G II'^{^),p G i7^~^(n) solving the problem 

'—/S.u + Vp = f inQ 

dYvu = h in Q (2 31) 

{pi - Bi{u))e3 = il; on^ ^ ' ' 

= on Sfe. 

Moreover, 

IMnr + Iblllf'-l ^ \\f\\Hr-2 + \\h\\Hr-l + ||V'|llfr-3/2 • (2.32) 

Proof See [5]. □ 

Lemma 2.16. Let r > 2. Suppose that f G H''-'^{Q),h G H''-^{Q),^p G H''-^I'^{Y:) and that 
h, (p satisfy the compatibility condition 



h= / (fi-es. (2.33) 

Then there exists unique u G II'^{Q,),p G II'^''^{0,) (up to constants) solving the problem 

' —Au + Vp = f inil 

div u = h in ^ , 

(2.34) 

u = if on S 

■u = on Sft. 



Moreover, 

Proof See [17]. □ 



^lll^^ + W^Pfur-^ < \\f\\lr-2 + \\h\\lr-l + Ij^'lj^.-V^ • (2-35) 
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2.3.2. The A— Stokes problems. In our local well-posedness argument, we will need the regular- 
ity of the solutions to the following ^—Stokes problem 

div j{ Syx{p,u) = in ft 

divA u = F^ inn 

Sa{p, u)M = F3 on E ^ ' ' 

u = on Sft. 

We shall use the regularity of (2.36) for each fixed t in the context of the time-dependent 
problem, and hence we will temporarily ignore the time dependence of 77, A, Af, etc, to view 
(2.36) as a stationary problem. 

Proposition 2.17. Let k > 4 be an integer and suppose that i] G H^'^^^'^iTi). There exists a 
universal £0 > so that if \\il\\k-i/2 — ^0, then solutions to (2.36) satisfy 

\Mr + Ibll.-i < + \\F\_^ + ||Fi,_3/2 (2.37) 

for r = 2, . . . ,k, whenever the right side is finite. 
In the case r = k + 1, solutions to (2.36) satisfy 

\Mk+i+\\p\\k^n-,+n+n_,/, 

/ \ (2-38) 

+ ll^ll.+l/2 (11^12 + 11^13+11^15/2)- 

Proof. Under the smallness assumption of r/, we can solve (2.36) by viewing it as a perturbation 
of the classical Stokes problem (2.31). See Proposition 3.7 of [12] for the detailed proof. □ 

2.3.3. The A—Poisson problem. We now consider the scalar elliptic problem 

'AAP = f^ inn 

p = f^ on E (2.39) 

yAP-v = f^ on Eb, 

where v is the outward-pointing normal on E(,. 

We first consider the weak formulation of (2.39). For this, we define a scalar T^} in a natural 
way through the norm 



^1= / J^Af?- (2.40) 



As in Lemma 2.9, this scalar norm generates the same topology as the usual scalar norm. 

For the weak formulation we suppose G {^H^{n)Y, f^ £ ^^/^(E), and G /f-V2(5]^). 
Let p G H^{n) be an extension of so that p = near E5. We switch the unknown to g = p — p 
and then define a weak formulation of (2.39) by finding a q £ ^H^{n) so that 

{q,f)w=-{P,^)w-{f\v), + {f,f)_y2 ioTall<fe'H\n), (2.41) 

where (•,-)^ is the dual pairing with ^H^{n) and (t)„i/2 dual pairing with //^/^(E^). 

The existence and uniqueness of a solution to (2.41) follow from standard arguments, and the 
resulting p = q + p € H^{n) satisfies 

\\p\\ni ^ ||/1(o//i(n))* + l|/^ll//i/2(s) + l|/^llH-i/2(Ef,) • (2-^^) 

In the event of determining the initial pressure the action of G {^H^(n))* will be given in 
a more specific fashion by 

{f\<p), = (90, v)no + (C, V^(^)^o (2.43) 
for go,G G H^{n) with H^ollo + IIG'llo = l|/^ll(Oi?i(Q))*- Then (2.41) may be rewritten as 

{VaP + G, Vav)ho = - {90, fho + {f, <^>_i/2 for all if G °iJ^(n). (2.44) 
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(2.45) 



\r-3/2 



This serves as motivation for us to say that p is a weak solution to the PDE 

Now we record the analogue of Proposition 2.17 for the problem (2.39). 

Proposition 2.18. Let k > 4 be an integer and suppose that rj G H'^^^/'^{T.). There exists a 
universal £q> Q so that if ||??||fc_]^/2 — ^O; then solutions to (2.39) satisfy 

lbll.<||/^L_2 + 11/11.-1/2 + 11/' 
for r = 2, . . . ,k, whenever the right side is finite. 

Proof. The proof is similar as Proposition 2.17. 

2.4. Transport estimate. Consider the equation 

{dtri + u ■ V^rj = g in S 
ri \t=o= m- 

We have the following estimate of the transport of regularity for solutions to (2.47). 

Lemma 2.19. Let rj be a solution to (2.47). Then there is a universal constant C > so that 
for any < s < 2 

sup Mr)\\Hs < exp (c [ ||V*n(r)||^3/2 dr^ (\\vo\\hs + Mr)\\Hsdr^ . (2.48) 



(2.46) 



□ 



(2.47) 



0<r<t 



Proof. Use p 

^3/2 



Bi^nL^ 





P2 = 2, N = 2, and a = s in Proposition 2.1 of [11] along with the embedding 

□ 



3. Local well-posedness of the k-problem 



The aim of this section is to solve the following regularized K-problem, where k > is the 
artificial viscosity coefficient, 



dfU — dtrjbKdsu + u ■ V^n — A^u + = in 



in Q 
on S 
on E 
on S(, 



(3.1) 



div_4 n = 

Sa{p, u)M = rjM- aHM 
dtrj — kA^t] = + u- M 
u = 

Xu,ri) |t=o= (-"c^/o)- 

Here ^ is the compensator determined by the initial data, and we postpone its construction in 
Section 3.3.2. We will temporally ignore the value of a since it plays no role in this section. 

3.1. The linear yl— Navier-Stokes problem. In this subsection we will consider the following 
linear problem for {u,p), where t] (and hence A,M, etc), and F'^ are given, 



(3.2) 



Since we are considering the same problem as Section 4 of [12], most of the arguments will 
be almost same as there. However, we will relax the assumptions on rj; as we will see, such 
relaxation will be crucial for our local- wellposedness of the K-problem (3.1) (and the original 
problem (1.11)) with the "semi-small" initial data. For the sake of completeness, we will repeat 



' dtu- Aau + VaP = 


in 


div^ u = 


in 


< SA{p,u)J\r = F^ 


on S 


u = 


on S;, 


U \t=0= Uq. 
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main arguments from [12]; along the context we will point out the differences and the necessary 
modifications due to our weaker assumption on ij. 

Although our aim is to construct strong solutions to (3.2) with high regularity, we need the 
definition of a weak solution to (3.2). Assuming the existence of a smooth solution to (3.2) and 
taking the L^L^ inner product with Jv for v G T-L^, we obtain 

{dtU,v)y^o^ + ^ {u,v)y^i^ - {p,diYAv)no, = {F\v)^o^ - (i^^^')o,s,T' (^•2) 

where (^OosT L'^L^ inner product in S x (0,r). If we were to restrict test functions 

to V E Xt (defined by (2.21)), then the term (]?, div^ t;)^o would vanish above, and we would 
be left with a "pressureless" formulation of the problem involving only the velocity field. This 
leads us to define a weak formulation without the pressure as follows. Suppose that 

F e {XtT and uq e y{0), (3.4) 

where ^(0) is defined by (2.19). Then our definition of a weak solution requires that u satisfies 

i {dtu, + i {u, = {F, for every G A't, (3.5) 
[«(0) = no, 

where (•, is the dual pairing with Xt- Note that the third condition in (3.5) makes sense in 
light of Lemma 2.12. Our formulation requires only that u G Ay, which means that F G {Xt)* 
is natural. Within the context of problem (3.2), F is most naturally of the form appearing on 
the right side of (3.3); if F admits a representation of this form, we may say that a solution to 
(3.5) is a weak solution of (3.2). We will not need to directly construct weak solutions to (3.5). 
Rather, weak solutions will arise as a byproduct of our construction of strong solutions to (3.2). 

Now we turn to the construction of strong solutions to (3.2). We will first construct the 
strong solutions with lower regularity by assuming that 

G ( [0, T] ; 1 (J^ ) ) n C° ( [0, T] ; (J^ ) ) , 

F=^GL2([0,T];i7=^/2(S))nC°([0,T];i7V2(s)), (3.6) 

dt{F^ - F^) G {XtT, and uq G H'^{n) n A'(O). 

The solution that we construct will satisfy (3.2) in the strong sense, but we will also show 
that Dtu satisfies an equation of the form (3.2) in the weak sense of (3.5). Here we define 

Dtu ■- dtu - Ru for R := diMM'^ (3.7) 

with M the matrix defined by (2.27). We employ the operator Dt because it preserves the 
div^ —free condition. Before turning to the result, we define the quantities 

JC{ri) := sup ( ||r?||9/2 + \\dtri\\y2 + ||^t^'?||5/2 ) '^(^) '= ^7/2 ■ (3-8) 

0<t<T ^ ' I " MO/^y o<t<T ' 

Theorem 3.1. Suppose that F^ ,F^ ,uq satisfy (3.6), and thatuo, F^{Q) satisfy the compatibility 
condition 

Uo{F^iO)+BAom^ro)=0, (3.9) 

where Uq is the orthogonal projection onto the tangent space of the surface {xs = rjo} defined 
according to 

Uov = v-{v J\fo)Mo \Mo\~^ . (3.10) 

Further suppose that /C(ry) < oo, and that C{r]) is less than the smaller of Eq from Lemma 2.9 
and Proposition 2.17. Then there exists a unique strong solution {u,p) to (3.2) so that 

ue XTnC{[0,T];H'^{n))nL^{[0,T];H^{n)), 

dtu e C{[0,T];H°{n)) n L'^{[0,T];H^{n)),Dtu e XT,d^u e (Xt)*, (3.11) 
peCi[0,T];H\n))nL\[0,T];H\^})). 
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The solution satisfies the estimate 

< P(l + /C(r?)) exp (P(l + /C(7?))T) (||uo||^ + \\F\0)\\1 + \\F^O)\\l^^ (3.12) 

The initial pressure, p{0) G H^{Cl), is determined in terms of uq, F^{0), F^{0) as the weak 
solution to 

'div^„(VA,p(0) - Fi(0)) = -diYAoiRoUo) G 

p{0) = (F3(0) +D^,KoAAo) • AAo lATop' G ijV2(s) (3.13) 
(V^oP(O) - F^O)) ■ ly = A^ouo ■ G iJ-V2(Sb) 
m i/ie sense of (2.45). Also, Dtu{0) = dtu{0) — RqUq satisfies 

Dtu{Q) = AaoUo - V^oP(O) + F\0) - Rouo G yiO). (3.14) 
Moreover, Dtu satisfies 

'dt{Dtu) - Aj^iDtu) + VAidtp) = DtF^ + m ft 

div AiDtu) =0 inn 

SAidtP, Dtu)M = dtF^ + on^ ^ ' ' 

DfU = on S;,, 

in the weak sense of (3.5), where G^,G^ are defined by 

G^ =-{R + dtJK)AAU - dtRu + {dtJK + R + R^)VaP 

+ dwA{^A{Ru) - R^AU + ^dtAu), (3.16) 

=nA{Ru)-^ - {pl - ^Au)dtM + Ba^^uA/". 

Here the inclusions (3.11) guarantee that G^ and G^ satisfy the same inclusions as F^,F^ listed 
in (3.6), whereas (3.13) guarantees that the initial data Dtu{Q) G 3^(0). 

Proof. The theorem is almost same as Theorem 4.3 of [12], and we shall first give a brief sketch 
of the proof and then explain the differences and the necessary modifications in our situation. 

The result will be established by first solving the pressureless problem in the weak sense of 
(3.5) via a time-dependent Galerkin method where the desired countable basis of Xt may be 
constructed with aid of Proposition 2.13. We then use Proposition 2.14 to introduce the pressure 
as a Lagrange multiplier, which gives a weak solution to (3.2). We then apply Proposition 2.17 
to improve the regularity of the solution so that the weak solution is indeed a strong solution. 

Along the proof the only difference in the assumptions of the two theorems is that Theorem 
4.3 of [12] requires that IC{r]) is sufficiently small, and our theorem only requires that part of 
IC{r]), C{r]), is sufficiently small but IC{rj) can be arbitrarily large. Note that the smallness 
of C{ri) is sufficient to guarantee the applications of the lemmas and propositions used in the 
proof of Theorem 4.3 of [12]. So the proof of Theorem 4.3 of [12] works in our present case, 
and the only difference resulting in the conclusions of the two theorems is that the polynomial 
P(l + IC{r])) appears in (3.12); in Theorem 4.3 of [12] since JC{rj) is small, such polynomial is 
replaced there accordingly by 1 + }C{ri) as P(l + }C{ri)) < 1 + /C(7?). □ 

Now we investigate the higher regularity of the strong solution obtained in Theorem 3.1. In 
order to define the forcing terms and initial data for the problem that results from temporally 
differentiating (3.2) several times, we define some mappings. Given F^,F^,v,q we define the 
mappings for forcing terms 

(3\v, q) = -{R + dtJK)AAV - dtRv + (dtJK + R + R^)VAq 

+ div a{^a{Rv) - R'S^AV + ^dtAv) on Q, (3.17) 

&^{v, q) = BA{Rv)J\f - {ql - BAv)dtJ^ + Bq^avM on E, 
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and the mappings for the initial data 

&^{F^,v,q) = Ajs,v -VAq + - Rv on 

f{F\.v) = {F^+'BAvAf)-J\f\J\f\~^ onS, 

fiF\v) = {F^ + Aav)-u onSb. 



(3.18) 



In the definitions of & and f we assume that A,M,R, etc are evaluated at the same t as 
F^,F^, V, q. These mappings allow us to define the forcing terms as follows. We write = F^ 
and = F^, and then we recursively define that for _7 = 1, . . . , 2N, 

i-i 

F^'^ := DtF^'^-' + <S\Dl-\,di-^p) = DlF^ + ^ ^^^©^(Dr^"'"^ ^"^"V) on Q, 

(3.19) 

e=o 

These are the forcing terms that appear when we apply j temporal derivatives to (3.2) (see 
(3.49)). In order to estimate these forcing terms, we define 

2N-1 „ 2N 

S{F\F^) := ^ \\diF' 

j=0 

2N-1 

2N-1 
j=0 



j=0 



l2fjiN-2j-l/2 



];^oafj4N-2j-2 



+ 



4N-2j-2 



+ 



dlF'iO) 



loafj4N-2j-3/2 
2 



(3.20) 



4N-2j-3/2 



Lemma 2.6 implies that if ^{F^,F^) < oo, then for j = 0, . . . , 2N — 1, 

diF^ e C{[0,T];H'^^-^^-'^{n)) and fi^F^ ^ ^([0, T]; iJ^^-2i-3/2(s)). (3.21) 

The same lemma also implies that the sum of the L°°H norms in the definition of d{F^,F^) 
can be bounded by ^o{F^ , F^) plus the sum of the L^H^ norms. For rj we define 

2N+1 „ 27V „ 



^iv) ■= \\v\\l2H4N+l/2 + E W^tV 



j^2jj4.N-2j+3/2 



+ hWlo.H^N +E 



J^<x,Jj4:N-2j + l/2 



2N 



4N-2j+l/2 



(3.22) 



Again, if ^{r]) < oo, then rj G C{[0,T]; H"^^ {^)) and dfrj G C([0, T]; ij4^-2j+V2(s)) for j 
1,...,2N. We also define 



2{ri) := sup ||r7||^4iv-i/2 and iio('7) := ||r/o||4Ar-i/2 • 

0<t<T 



(3.23) 



Note that JC{r]) < ^iv) C{r]) < £,{r]), where IC{r]) and C{r]) are defined by (3.8). 

We now record an estimate of F^'^ in terms of ^{F^,F^), ^{rj) and certain norms of u,p. 



Lemma 3.2. For m = 1, . . . , 2N — 1 and j = 1, . . . ,m, 

\\F''ml2H2m-2J+, + \\F^''\\l2H2m-2J+S/2 < P(l + 5(ry)) (^5(^1 , F^) 



i-1 



+EF 



u 



1,2 ff2m-2j+3 



+ 



dtu 



l^ao u2m-2j+2 



+ 



l2fj2m-2i+2 



+ 



dip 



(3.24) 
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I \\L°°H'^m-2j 



u 



2 






+ 


dip 




J^oo Jj2m — 2j + 2 


l^ooJ{2m-2j + l ; 



(3.25) 



and 



\dtiF 



l.m 



7l3,m^ 



+ 



m— 1 

E 

£=0 



+ 



9^ 



2 



+ 



dip 



+ 



dip 



(3.26) 



Similarly, for j 
\\F''^{0)\\ 



2 

4Ar- 



2j-2 



+ 



2iV- 1, 

F''^m 



2 

4A?"- 



■2j-3/2 
J-1 ., 



< P(l +;?o(r?)) ^o(i^\i^=^) + ^ pfn(O) 



£=0 



2 

4:N-2e 



+ 



dfpiO) 



2 

4Af-2£- 



(3.27) 



Proof. The lemma follows from simple but lengthy computations, invoking standard arguments. 
Although our lemma is similar as Lemma 4.5 of [12], we need to point out the differences and 
present the necessary modifications of the proof there. The first point is that we assume ^{rj) 
to be finite but not small, so the polynomial P appears in our estimates. The second point is 
that some norms in our 5^(t/) are weaker than those of ^(rj) in Lemma 4.5 of [12], but a careful 
check that we shall explain below would reveal that they are sufficient to guarantee the validity 
of these estimates. Hence, our estimate is a refined version of Lemma 4.5 in [12] in some sense. 

We use the definition of F^'^ , F^'^ given by (3.19) and expand all terms using the Leibniz rule 
and the definition of Dt (defined by (3.7)) to rewrite i^*'-^ as a sum of products of two terms: 
one involving products of various derivatives of fj, and one linear in derivatives of u, p, F^, or 
F^. We only focus on the fj part to see that it can be bounded by our polynomial, while the 
other part follows in the same way as Lemma 4.5 of [12]. Note that the highest derivatives of f/ 
appearing in F^'^ are div^fj and di~^^Vf], while the one appearing in F^'^ is c^V^f/. To prove 
(3.24), by Lemma 2.1 and the trace embeddings, we have that for j = 1, . . . , 2N — 1, 



< 



2 



23-1 



< 



L2H4N- 
2 



23+3/2 



]^2ff4N-2j 



< 

-1/2 (S) ~ 



L2H4N- 
2. 2 



2J-1/2 



L2ij4JV-2j 



L2H*N-2(j + l)+3/2 



<d{v), and (3.28) 



^2^4JV-23+3/2 



With (3.28) in hand, we may conclude (3.24) as in Lemma 4.5 of [12] by estimating the re- 
sulting products by using Lemma 2.2 in conjunction with the usual Sobolev embeddings, trace 
embeddings and Lemma 2.1. The estimates (3.25) and (3.27) follow from similar arguments. 

To prove (3.26) , we note that the highest derivatives of rj are d^^V^fj and d^^'^^Vfj appearing 
in dtF^'"^ and 9|^V^f/ appearing in dfF^'"^ when m = 2N — 1. We remark that the arguments 
for (4.78) in Lemma 4.5 of [12] or the way of (3.28) can not be applied to estimate (3.26) due to 
our weaker d{v)- We must resort to use the structure of the estimated terms. First, to control 
d^^'^'^Vfj in n, we have that for all ip G o^^(^^), 



Jn 



(3.29) 



Then by the duality, the trace embedding and Lemma 2.1, we obtain 



L2(oHi(n))" 



< 



L2//-1/2 



+ 



< 



L2H0 



L2//-1/2 



(3.30) 
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On the other hand, we have to treat together the terms in the dt{F^''^^~^ — i?3,2Af-i^ that 
leads to df^V^f] in Q and d'^^V'^fj on S. In view of (3.19) for j = 2N — 1, these terms 
are div^(B^(9^^~^i?M)) in 0, and —I])ji^{d^^^^Ru)J\f on S, respectively. We need to estimate 



Then by the duality wc have 



. Indeed, for il) G Xt, we observe that 



dwA{^A{dr~^Ru)) - ^A{dr-'Ru)N 



2N-1 



< 



Ru 



To estimate the right hand side of (3.32), we may use the bound, by Lemma 2.1, 

Ii5rv^iil.^i<ii9r^iil.^3/.<;?(r/). 

Then with (3.30), (3.32) and (3.33) in hand, we can conclude (3.26) as that for (3.24). 

Next we record an estimates for the difference between dtv and DtV for a general v. 

Lemma 3.3. If k = 0, . . . , AN — 1 and v is sufficiently regular, then 

\\dtv - DtvWl^H, < P{1 + m) IMhm , 
and ifk = 0,...,4:N-2, then 

// m = 1, . . . , 2A'' — 1, j = 1, . . . , m, and v is sufficiently regular, then 



&>v-Dl 



(\\dfv 


2 

+ 











]^oo [j2m-2j + 2 



d^v-Di 



Loofj2m-2j+2 



Loofj2m-2j+2 



and 



<p(i+ 5(^7)) Elk* 



+ 



im+l„ 



Also, if j = 0, . . . , 2N, and v is sufficiently regular, then 

div{0) - DlviO) <P{l+Uv))J2\rtv{0) 

4N-2J ^ II 



2 

AN-2e ' 



(3.31) 
(3.32) 

(3.33) 
□ 

(3.34) 
(3.35) 

(3.36) 

(3.37) 



(3.38) 



(3.39) 

□ 



Proof. The proof is similar to that of Lemma 3.2, and is thus omitted. 

Lemma 3.4. Suppose that v, q, G^, are evaluated att = and are sufficiently regular. For 
j = 0, . . . , 2N — \, we have 



\\e\G\v,, 

<pii + \\vmlN + \\dtvm\lN-3/2){\\ + ll^i4Ar-2,-2) • 

If j = 0,...,2N -2, then 



(3.40) 



<p{^ + \\vmlN + \\dtvmlN- 



-3/2) 



.1||2 

l47V-2j-2 



•,3||2 

l4Ar-2j-3/2 



I' 

\4N-2j) ■ 

(3.41) 
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For j = 2N — 1, if div^(o) v = ^ in ^, then 

\\f{G\v)\\l^ + \\f{G\v)\\'_^^^ < P{1 + WmWrn) {\\gX + \\G%, + \\v\\l) . (3.42) 
Proof. The proof is similar to that of Lemma 3.2. We note that to derive the estimate of 

II 1 1 2 II 1 1 2 

(3.42), we may use the bound || A_4(o)f • ^\\h-i/2(^j^^-j ^ ||^-4(o)'"||o provided by Lemma 2.11 since 
div_4(o) ^A{o)''^ = 0) which is imphed by that div^(o) ^ = 0. □ 

Now we assume that uq G H^^{n), rjo G H'^^{T,), ^o{F^,F^) < oo and do{r]) < oo, and that 
S-oiv) ^ 1 (as defined in (3.23)) is sufficiently small for the hypothesis of Proposition 2.18 to 

1 1 2 1 1 2 

hold when k = AN . Note, though, that we do not need ||?/o||4jv ^^id IpolUiv to be small. We will 
iteratively construct the initial data D-j.u{0) for j = 0, . . . , 2N and dfp{0) for j = 0, . . . , 2N — 1. 
To do so, we will first construct all but the highest order data, and then we will state some 
compatibility conditions for the data. These are necessary to construct D^^u{0) and 5|^~^p(0), 
and to construct high-regularity solutions in Theorem 3.5. 

We now turn to the construction of Dlu{0) for j = 0, . . . , 2A^— 1 and c^p(O) for j = 0, . . . , 2N— 
2. For J = we write ^^'^(O) = ^^(0) G H^^-'^, F^'^{0) = F^{0) G ^^^^ D^u{0) = 

no G H^"". Suppose now that F^'^ G /f4iV-2£-2^ p3,e ^ ^4iV-2^-3/2^ j^e^^Q^ ^ jj4N-2e 
given for < £ < j G [0,2A^ - 2]; we will define d{p{0) G ij4Af-2i-i ^^^j d{^^u{0) G 

jjm-2j-2^ ^l,j+l(0) G i/4Af-2i-4^ F3'i+1(0) G H^N-2j-7/2^ ^j^j^j^ ^gg^g 

said data via iteration. By virtue of estimate (3.41), we know that = f^(F^'-'(0), £)^tt(0)) G 

H4N-2j-3^ j2 ^ f^F^-'iO),Diu{0)) G ii-4iV-2,-3/2^ f3 ^ f (^1^(0), Z?^n(0)) G ii-4iV-2,--5/2_ 

This allows us to define ^p(O) as the solution to (2.39) with this choice of f^,f^,f^, and then 
Proposition 2.18 with A: = 47V and r = 4iV - 2j - 1 < /c implies that ^p(O) G ij4iv-2j-i_ -^^^ 
the estimates (3.27), (3.39), and (3.40) allow us to define 

Di+\{0) := 0'{F^'^{O),DiuiO),d{piO)) G H^^'^^-^ 

^iJ+i(0) := DtF^'^iO) + <S\Diu{0), d{p{0)) G h'^^-^^'^ and (3.43) 
F^'j+\0) := dtF^'^iO) + (3^{Diu{0),d{p{0)) G ii-4iV-2j-7/2_ 

Hence, we have iteratively constructed all of the desired data except for D^^u{0) and d^^^^p{0). 

By construction, the initial data L>^ti(0) and dlp{0) are determined in terms of uq as well as 
d^F^{0) and d^F^{0) for ^ = 0, . . . , 2A^-1. In order to use these in Theorem 3.1 and to construct 
Df^u{0) and d^^~^p{0), we must enforce compatibility conditions for j = 0, . . . ,2N — 1. For 
such j, we say that the j*'* compatibility condition is satisfied if 

fL>M0)e-^(0)niJ2(Q) 

\no(F3.J(0) + D^D^n(O)Aro) = 0. ^ ' ' 

Note that the construction of Dlu{Q) and dlp{0) ensures that -D|m(0) G H^{Vt) and that 
diYAo{Dlu{Q)) = 0, so the condition L'^n(O) G Af(0) n H'^{Q.) may be reduced to the condi- 
tion r»|n(0)|s, = 0. 

It remains only to define ^^^^"^(0) G and D^^u{0) G H^. According to the j = 27V - 1 
compatibility condition (3.44), div^g D^'^~^u{0) = 0, which means that we can use the estimate 
(3.42) to see that /2 = fiF^'^^-^iO),D^''-\{0)) G H^^ and = f (f1'2^-1(0), ^^^-^^(O)) G 
H-^/^. We also see from (3.44) that if we define 50 = - div^o(i?oA^^"^^i(0)) then go G 
Then, owing to the fact that G = — G H^, we can define d^'^~^p{0) G as a weak 
solution to (2.39) in the sense of (2.45). Then by the estimate (3.40) we define 

L>2^n(0) = ©0(Fi'2^-i(0), L>2iV-i^(o), a2^-ip(0)) G iJ^ (3.45) 
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In fact, the construction of df p{0) guarantees that Df u{0) € 3^(0). In addition to providing 
the above inclusions, the bounds (3.27), (3.41), (3.40) also imply the estimate 

2Af-l 



27V 
j=0 



4N-2j 



j=0 



+ ^ tep(O) ' <Pil + Uri))(\\uo\\lN+UF\F^)). (3.46) 

. — i II 4Ar— 2j— 1 V / 



Note that, owing to (3.39), (3.46) also holds with ^'u(O) replacing Dju(O) on the left. 
To state our result on higher regularity of solutions to the problem (3.2), we define 



2N 



2N-1 



Idju 



3=0 



L2H4N-2j+l 



+ 



+ E pip 



3=0 



L2H4N-2j 



2N 
j=0 



u 



l^ooH4N-2j 



2N-1 

+ E 

j=0 



(3.47) 



\dip 



lj00ff4N-2j-l 



Theorem 3.5. Suppose that uq e H'^^{n), ryo G i^'^^(S), d{F^,F^) < oo, ^iv) < oo and that 
£(ry) (as defined in (3.23) ) is less than the smaller of eo from Theorem 3.1 and Proposition 2.18. 
Let Diu{0) e H'^^-'^i{Q) and d^^p{0) G ij4Ar-2j-i for j = 0, . . . ,2N - 1 along with Df^u{<S) G 
3^(0) all he determined as above in terms of uq and d'jF^iO), d^F^iO) for j = 0, . . . ,2N - 1. 
Suppose that for j = 0, . . . , 2N — I, the initial data satisfy the j^^ compatibility condition (3.44). 

There exists a T = T{^[r])) > so that if < T < T, then there exists a unique strong 
solution {u,p) to (3.2) on [0, T] so that 

d{u e C([0, T];H^^-^^{n)) n L\[0, T];H'^^-'^^+\n)) /or i = 0, . . . , 2N, df^+\ G {XtT , 
dip G C([0,T];if4^-2^-i(f])) n L\[0,T];H''^-'\n)) for j = 0, . . . ,2N - 1. 



(3.48) 



The pair {Dlu,dlp) satisfies the PDE 

(dt{D{u) - A_A{Diu) + VAidip) 
I div A{Diu) = 
\sAidip,Diu)M = F^'^ 
[Diu = 



p^,3 



in ri 
in 

on S 
on Sfo 



(3.49) 



in the strong sense with initial data {D^u{0) , p{0)) for j = 0, . . . ,2N — 1, and in the weak 
sense of (3.5) with initial data D^^u{0) G 3^(0) for j = 2N. Here the vectors F^'^ and F^'^ are 
as defined by (3.19). Moreover, the solution satisfies the estimate 

^{u,p) < P{l+Mv) + ^iv)) exp (P(l + d{v))T) {\\uo\\l^ + MF\F^) + 5(F\ F=^)) . (3.50) 

Proof. The proof is same to that of Theorem 4.8 of [12]. We will only provide a brief sketch of the 
idea of the proof. For full details we refer to [12]. The estimate (3.46) gives control of the initial 
data, the estimates (3.24)-(3.26) in Lemma 3.2 provide control of the forcing terms F^'^,F^'^ 
in terms of F^ , F^ ,r],u,p, and the estimates (3.34)-(3.38) in Lemma 3.3 provide control of the 
difference between dtv and DfV. Note that the smallness of £,{r]) is sufficient to guarantee the 
applications of the lemmas and propositions used in the proof of Theorem 4.8 of [12] so that the 
proof of Theorem 4.8 of [12] works in our present case. These and the j*'* compatibility condition 
(3.44) for j = 0, . . . , 2N—1 then allow us to iteratively apply Theorem 3.1 and Proposition 2.17 
to find that (D^n, ^p) solve (3.49) and satisfy the estimates (3.50). □ 

3.2. The linear regularized surface K-problem. In this subsection we will consider the 
following regularized surface problem with the artificial viscosity term and the compensator: 



V \t=o= Vo, 



in S 



(3.51) 
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where k > is the artificial viscosity parameter and the function ^ is the so-called compensator 
that will be determined below. We define 



2 

l,2H4,N 



2N 

5(f4) ■.= ^\\dlF 

j=0 

2N-1 



2N-1 

... + E W^."' 

j=0 



l^oofj4N-2j 



(3.52) 



4iV-2j-l 



Again, Lemma 2.6 implies that if ^^(F^) < oo, then g C{[0,T]; H^'^-^^'^J:)) for j = 

0, . . . , 2A'^ — 1. To state our result on the solutions to the problem (3.51), we define 



2Af+l 
j=0 



in 



]^2jj4N-2j+2 



2N 

i=o 



£ocJj4N-2j+l 



(3.53) 



Now we turn to the construction of the compensator function ^. In order not to break the 
compatibility conditions required for the initial data (no,?7o) of the nonlinear system (1.11) (i.e. 
(3.1) with K = 0), we need to require that the initial conditions for (3.51) are same as the one 
with K = 0. This is the reason that we introduce ^ here. Formally, we use the equation (3.51) 
with k = to construct the initial data (^r/(0), and then choose the appropriate function \1' 
so that the equation (3.51) with A; > takes the same initial conditions as (3.51) with A; = 0. 
We construct ^' as follows. We assume that r]o G i?^^+i(S), and we construct the initial data 
^r?(0) for i = 1, . . . , 2N as dlr]{0) := cP^F*{0). Now we define tj^ to be the extension produced 
by Lemma 2.8 that achieves the initial data (^r/(0) for j = 0, . . . , 2N. If we define ^' = — A^ry*^, 
then (?j*(0) = — A*c^?7(0). This guarantees that we essentially add nothing in (3.51) at the 
initial time t = 0. Moreover, we have 



2N 
j=0 



L2H4N-2j 



2N 
3=0 



2N 



< 



J^x^4N-2j-l 



3=0 



AN-2j+l 



(3.54) 



holl4JV+l +5^o(^^ 



Thanks to this compensator, the constructed data dlrf{0) for j = 1, . . . ,2N can serve as the 
initial data for (3.51) for all k. 

Theorem 3.6. Suppose that r]o G i?^^+i(S) and that ^r?(0) for j = 1, . . . , 2N are determined 
as above. Then the problem (3.51) admits a unique solution rj that achieves the initial data 
^7/(0) for j = 0, . . . , 2N and satisfies 



m < C« (||??o||4iv+i + MF") + d{F' 



(3.55) 



Proof. The solvability of (3.51) is classical, and we only focus on the derivation of the estimate 
(3.55). The standard energy estimates on (3.51) together with Cauchy's inequality yields 



d II ||2 , IIY7 ||2 ^ lhT(l|2 

^ ll^lUiv+i + II V *?7ll47v+i ~ '^II^IUat 
Integrating the above in time, by (3.54) we obtain 



K 



4||2 



(3.56) 



kllL°°H4JV+l + K \\r]\\1^2ff4N+2 < ||??(0)||4JV_,_i + K ||^'|||2^4JV H 1 1 1 1 L2^f4iV 

<\\rio\\lN+i+UF') + l^iF% 

rv 



(3.57) 
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On the other hand, we directly use the equation (3.51) together with (3.54) to have that for 
j = l,...,2Ar + l 



< 



K 



4 



]^2ff4N-2j+4 
2 

j^2ff4N-2U-l)+2 



l2fj4N-2j+2 



+ 



j^2fj4N-2j+2 



+ ll%ll47V+l 



A simple induction on (3.58), together with (3.57), yields 

2Af+l 



< 



l^2ff4N-2j+2 



\r]\\]^2jj4N+2 + II'VoIUjv+i 



(3.58) 



(3.59) 



<\\m\\lN+i + UF^) + ^{F^)- 

This and (3.57) yield (3.55) by further applying Lemma 2.6 to bound the sum of the L°°H^ 
norms in the definition of ^{rf) by Ht/oIUjv+i + i?o(^'^) plus the sum of the L^H^ norms. □ 

3.3. Preliminaries for the nonlinear /c-problem. We will eventually use our linear theories 
for the problems (3.2) and (3.51) in order to solve the nonlinear problem (3.1). To do so, we 
define forcing terms F^,F^ to be used in the linear theory of (3.2) and to be used in the 
linear theory of (3.51) that match the terms in (3.1). That is, given u,r], we define 

F^{u, T]) = dtfjbKdsu - u ■ V _au, F^{r)) = rjAf - aHAf, and F^(u, r]) = u ■ N . (3.60) 
where A,Af, K, H are determined as usual by rj. 

3.3.1. Data estimates. In the construction of the initial data performed in the previous two 
subsections, when constructing 5^u(0) for j = 1, . . . ,2N it was assumed that &lr]{0) for j = 
0, . . . , 2A'" and dlF^{0), c^F'^(O) for j = 0, . . . , 2N — 1 were all known; when constructing (^?7(0) 
for _7 = 1, . . . , 2N it was assumed that ^F^(O) for j = 0, . . . , 2N — 1 were known. Since for 
the full nonlinear problem the functions u,p,ri are unknown and their evolutions are coupled 
to each other through F^(u,ri), F'^[rj) and F'^{u,ri), we must revise the construction of the 
data to include this coupling, assuming only that uq and rjo are given. This will also reveal 
the compatibility conditions that must be satisfied by uq and rjo in order to solve the nonlinear 
problem (3.1) (and (1.11)). To this end we define 



So-.-- 



1 1 2 1 1 2 

\uoUn + I^oIUat+i 5 and £0 - 



|%|l4Ar-l/2 • 



For our estimates we must also introduce 

2Af 

(&o{u,p) ■.= ^\\diu{0) 

j=0 

and 



2A-1 



4Ar-2j 



2N 



j=0 



j=0 



4Af-2j+l 



4Af-2j-l 



(3.61) 



(3.62) 



(3.63) 



We will also need a more exact enumeration of the terms in ^q{u,p), ^o{rj), ^q{F^ , F^), and 
doiF^) (as defined in (3.20) and (3.52), respectively). For j = 0, . . . , 2Ar - 1 we define 



and for j 



^i{F\u,v),F\r,)) ■.= J2\\dfF\0) 

di{FHu,r,)) := \\diF\0) 
0, . . . , 2iV we define 



4iV-2^-2 
2 



+ 



dfF\0) 



iN-2i-Z/2 



4N-2j-l 



4Ar-2^+l 



(3.64) 
(3.65) 

(3.66) 
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and 



2 

4N-2e 



i-1 



4N-2e-l 



(3.67) 



I l|2 



Here in (3.67) for j = we mean £q(u,p) := 

We record the estimates of these quantities that are useful in deahng with the initial data. 



Lemma 3.7. For j 

For j 



0, ...,27V- 1, it holds that 

di{F\u,v),F^{7^)) < P,i<£l^\r]),<tiiu,p)). (3.68) 

= l,...,2iV - 1 let Fi'J(0) and F^'^{0) be determined by (3.19) and (3.60), using 
dfri{0), dfu{0), and dfp{0) for appropriate values of £. Then 



l^''^(0)|lL-2,-2 + ll^''^(0)|l4iV-2,-3/2 < Pji^l^\rj)^ ^ii^^P))- 



4N-2j 



<pj{4{v),4{u,p)). 



(3.69) 
(3.70) 
(3.71) 



Forj = 0,..., 2N it holds that 

diuiO)-Diu{0) 

Forj = Q,...,2N-l it holds that 

diiFHu,rj)) < P,iiBiirj),<&iiu,p)). 
Here Pj{-,-) o, polynomial so that Pj(0,0) = 0. 

Proof. These bounds may be derived by arguing as in the proof of Lemma 3.2. As such, we 
again omit the details. □ 

This lemma allows us to modify the construction presented in the previous two subsections to 
construct all of the initial data dl it(0), dl r]{0) for j = 0, . . . , 2N and &lp{0) for j = 0, . . . , 2N— 1. 
Along the way we will also derive estimates of <Eo{u,p) + G;o('?) in terms of £q and determine 
the compatibility conditions for uo,rio necessary for existence of solutions to (3.1) (and (1.11)). 

We assume that uo,r)o satisfy < oo and that £o < 1 is sufficiently small for the hypothesis 
of Proposition 2.18 to hold when k = AN. As before, we will iterativcly construct the initial 
data, but this time we will use the estimates in Lemma 3.7. Define dtrj{0) := F'^{0) = uq • A/q, 
where uo G H'^^^^^'^{T,) when traced onto S. Then ||9try(0)||4jY_;^y2 ^ "^0(1 + £o), smd hence in 
particular that <Eq{u,p) + ^^{ri) < P{£o). We may use this bound in (3.68) with j = to find 

^^{F\u,r,),F\ri)) < Po«(«,p), ^^(t?)) < P(^o). (3.72) 

Suppose now that j G [0, 2A^ — 2] and that dfu{0) are known for £ = 0, . . . ,j, dfri{0) are 
known for £ = 0, . . . ,j + 1, and dfp{0) are known for £ = 0, ... ,j — 1 (with the understanding 
that nothing is known of p{0) when j = 0) , and that 

4{u,p) + 4-^\ri) + di{F\u,rj),F^{u,ri)) < P{£^). 

According to the estimates (3.69) and (3.70), we then know that 



|F^'^(o)ii: 



+ 1 1 ^^'^'(0)1 



l47V-2j-2 ' ir '^"^ll4iV-2jr'-3/2 

By virtue of estimates (3.41) and (3.73), we know that 



+ 



4N-2j 



< P{£o). 



(3.73) 

(3.74) 



fHF^'^(0),i^^(0)) 



+ 



f{F^^\Q),Dlum 



AN-2j-3 
2 



+ 



fiF''^{0),Diu{0)) 

< P{£o). 



4N-2j-3/2 



(3.75) 



4Ar-2j-5/2 

This allows us to define &lp{0) as the solution to (2.39) with given by f-^, f^, f^. Then 

Proposition 2.18 with k = AN and r = AN — 2j — 1 < k implies that 

2 



dip{0) 



4N-2j-l 



(3.76) 
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Now the estimates (3.40), (3.73), and (3.74) allow us to define 
and owing to (3.70), we have the estimate 



Now we define 



4iV-2(j+l) 



i+1 



£=0 

The estimate (3.71) together with (3.73) and (3.78) then imply that 



< P{£o) 



4Af-2(i+2)+l 

We may combine (3.73) with (3.76)-(3.80) to deduce that 

but then (3.68) implies that '^■{^^{F^{u,ri),F^{u^rj)) < P{£q) as well, and we deduce that the 
bound (3.73) also holds with j replaced by j + 1. 
Using the above analysis, we may iterate from j 



(3.77) 
(3.78) 

(3.79) 

(3.80) 
(3.81) 



0, . . . , 2A'" - 2 to deduce that 



-i(n,p) + {rj) +r^'-\F\u,r^),F-\u,rj)) < P{£o). 



-2N-1 



(3.82) 



After this iteration, it remains only to define ~ p{0) and D^u{0). In order to do this, 
we must first impose the compatibility conditions on uq and rjo. These are the same as in 
(3.44), but because now the temporal derivatives of t] have been constructed as well, we restate 
them in a slightly different way Let ^w(O), F^'^(0), F^'^{0) for j = 0, . . . , 2iV - 1, d{r]{0) for 
j = 0,. . . , 2N, and 9/p(0) for j = 0, . . . , 2N — 2 be constructed in terms of uq as above. Let 
Ho be the projection defined in terms of tjo as in (3.10) and Dt be the operator defined by (3.7). 
We say that uo, rjo satisfy the (2iV)*'* order compatibility conditions if 



Dlu{0) = 



in Q 

on Ef, 



(3.83) 



_no(F3'J(0) + BaoDIu{0)Mo) = on S 

for j = 0, . . . , 2A'^— 1. Note that if UQ,r]Q satisfy (3.83), then the j*'^ order compatibility condition 
(3.44) is satisfied for j = 0, . . . , 2iV - 1. 

Now we define ^(^^-^^(0) and ^^^^(O). We use the compatibility conditions (3.83) and argue 
as above and in the derivation of (3.42) in Lemma 3.4 to estimate 



:2,rf,2iV-l(o)^^2iV-l^(Q)) 



f{F' 



1/2 



+ 



,3,^l,2iV-l(o)^^2iV-l^(Q)) 



f{F^ 



-1/2 



and 



nl,2Ar-l/ 



(o)||;;+ divAoiRoDr-'m) <p{£ 



(3.84) 



(3.85) 



We then define d^^ ^p(O) G as a weak solution to (2.39) in the sense of (2.45) with this 
choice of f = f, f = f, go = - dwAo{RoDt^'^u{0)), and G = -F^'^^-^{0). The estimate 
(2.42), when combined with (3.84)-(3.85), allows us to deduce that 

32W-1 



(3.86) 



Then by (3.40) we set D^^u{0) = e^{F^'^^-^{0), D^^-\{0),d^^~^p{0)) G H^. In fact, the 
construction of df^~^p{0) guarantees that Df^u{0) G 3^(0). Arguing as before, we also have 

\\d^''uiO)\\l<P{£o). (3.87) 

This completes the construction of the initial data, but we will record a form of the estimates 
(3.82), (3.86)-(3.87) in the following proposition. 
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Proposition 3.8. Suppose that no,% satisfy £q < oo and that £,q < Eq is sufficiently small for 
the hypothesis of Proposition 2.18 to hold when k = AN. Let d^u{0), d^r]{0) for j = 0, . . . , 2N 
and dlpifS) for j = 0, . . . , 2N — 1 be given as above. Then 



£o<^oiu,p) + Mv)<P{£Q)- 

Proof. The proposition directly follows by summing (3.82) and (3.86)-(3.87). 



(3.88) 
□ 



Remark 3.9. Along the arguments above, it is easy to see that rjo £ iJ^^+-'^/^(S) replacing 

r]o G H^^'^^{I1) is still sufficient to guarantee the construction and that the constructed initial 
data can enjoy the following improvement, one may also refer to Section 5.2 of [12], 



27V+1 



(Bo{u,p) + \\vm\lN+i/2+ E IftViO) 



4iV-2j+3/2 



< P (ll^o|l4Ar + lkoll4Af+l/2 



(3.89) 



Note that we have also included the case j = 2N + 1 in (3.89). Indeed, using the geometric fact: 
dt{u ■ J\f) = DfU ■ M onTi, by Lemma 2.11 and since D^^u{0) € 3^(0), we have 

2 



_l/2 " ^ — ^011^-1/2(2) 

< P (ll«o||4Ar + ll%ll4Ar+l/2) • 



< D., 



2N„ 



^n(0)|| 



HO 



(3.90) 



3.3.2. Construction of compensator function. Let dfr]{0) for j = 0, . . . , 2N be determined as 
in Section 3.3.1. Now we define ^ to be the extension produced by Lemma 2.8 that achieves 
the initial data — A^=(^?7(0) for j = 0,. . . ,2N. Then we essentially add nothing on the fourth 
equation of (3.1) at the initial time. Moreover, by Lemma 2.8 and (3.88), we have 



2N 
j=0 



2N 



l2H4N-2j 



j=0 



J^oofj4N-2j 



2N 

<E||«o) 

j=0 



4N-2j+l 



= ^oiv) < P{£o). (3.91) 



Thanks to this compensator, the initial data {uq, 779) a-nd those initial data actually constructed 
for the problem (1.11) can serve as the initial data for the problem (3.1). 

3.3.3. Forcing estimates. We will now estimate various norms of Fi(n,77), F^{r]) and F^{u,rj) 
that are contained in ^{F^,F^) and d{F^) (as defined in (3.20) and (3.52), respectively) in terms 
of the certain norms of u and rj. It is a very key point to notice that the rj and u appearing in 
^{rj) and ^{u,p) enjoy the higher regularities than that required in estimating the norms of the 
F^ terms. To state the estimates, we define 



2N 

d2N{u)=^\\dl, 
j=0 



]^2fj4N-2j+3/4 



2N-1 

+ E IH'^ 

j=0 



];^oofj4N-2j-l/4 



(3.92) 



and 



2N 
j=0 



L2H4N-2j+l 



LoaJj4N-2j 



(3.93) 



2N-1 

. + E Ih' 

j=0 

We shall recall the definitions of A2n{u), .^2iv(?7), ^{u,p), div): ^'^d A{r]) from (2.12), (2.14), 
(3.22), (3.47) and (3.53), respectively. Note that 5^2Ar('") < ^2n{u) < A{u,0) and d2N{v) < 
^2Niv)<^{v),Siv)<^iv)- 



Lemma 3.10. Suppose that ^{rj) < 00 and S^2Ar('u) < 00. Then 
d{F\u,ri),F'{v)) < P{l + m) im + 

and 

^{F'^iu,!!)) < P{l+^2N{ri))d2Niu). 



(3.94) 
(3.95) 
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Proof. The proof is similar to that of Lemma 3.2. We only estimate the terms IL2(oifi(n))* 

and ||-F^||^2^4Jv-i/2 in the definition of 'S{F^,F^) since they are a bit more involved and need 
more care, and the other estimates follow similarly and simpler. 

We note that the highest derivatives appearing in df^F\u,r]) are d^^+^fj and d'^^Vu. The 
q2N+i- j^g^g been already controlled in (3.30). For the df^Vu, we observe that for (p G qH^{^1), 

{d^Vu,^)^ = (9r-,V^)o,E- (^r-,Vv.)o,^. (3.96) 

Here (•, •)^ is the dual paring in qH^{P), (•, •)q and (•, •)q is the LP' inner product in S and il, 

respectively By the duality and the trace embedding that H^/^{Q) ^ H'^/'^{T) C ii"~^/2(S), 
we obtain 

W^4lW^)Y - W^t'^A^L^H^/^ + W^L^HO < ||5r^||'2^3/4 < ^2n{u). (3.97) 
With (3.30) and (3.97) in hand, we may conclude as in Lemma 3.2 that 

<p{i+m){m + {^2N{v)f). 

On the other hand, note that the highest derivative appearing in F'''{ri) is V^??, and we have 

\\^lv\\\2HiN-l/2 < ||?7lli2^4;V+3/2 < ^{t]) ■ (3.99) 

We may then conclude as in Lemma 3.2 that 

||i^i'2H4iv-i/2 < ^(1 + ^{r])Mr])- (3.100) 
Consequently, estimating the other terms in the similar way, we deduce (3.94)-(3.95). □ 

Remark 3.11. Note that at this moment there is no power of T in front of the estimates 
(3.94) -(3.95). The key point is that in the iteration argument in the next subsection we can get 
better estimates of rj and u, say, ^{rj) and A{u, p) . Eventually, this will make the appearance of 
a power of T in front of some estimates in the iteration argument, which allows us to achieve 
a local well-posedness without requiring all the norms of initial data he small. 

3.4. Local well-posedness of the nonlinear K-problem. In order to construct the solution 
to (3.1), we will pass to the limit in a sequence of approximate solutions. In this subsection we 
will explicitly enumerate the various generic polynomials appearing in estimates as Pi, ... , Pio 
so that they can be referred, but they will be allowed to depend on k. Before proceeding further, 
we recall from the previous subsection that 

€o{u,p) + <to{ri)+^o{F\u,ri),F^{ri))+:So{F\u,7i))<Pi{£o). (3.101) 

We now construct the sequence of approximate solutions as follows. First, we extend the 
initial data df.u{0) G H^^^'^^Q) to the time-dependent functions so that dlu^{0) = dlu{0) 
and extend the initial data ^^//(O) G //^^^^J'+^(S) to the time-dependent functions r]^ so that 
djrj^{0) = dlr]{0). We do so by applying lemmas 2.7-2.8 respectively, and we may combine the 
estimates of A2n{'u'^) and ^2n{v^) provided by lemmas 2.7-2.8 with (3.101) to bound 

J^2iv(^°) + i^27v(??°) < ^2(^0). (3.102) 

This is our start point for constructing {u"^ , p"^ , rj™-) for m > 1. 

Now assuming that (ti™,??'") are known and satisfy certain estimates, we then construct 
^^m+i^pm+i^^m+i^ as the solution to 

'^j^m+i _ A_4™+in™+i + V^m+ip™+i ^ Fi(tx"», j7"^+i) in n 

div_4m+i = inn 

Sa-- (p™+^ n™+i)AA'"+i = F3(^m+i) J. 

5^^m+i _ = + F^{u"\ 7?"^) on E 

u"'+^ = on 



(3.103) 
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Note that our construction is in its special way; we first use (m™,?]™) to construct r/™"*"^, and 

We shall derive the uniform estimates for 



then we use (m™,?7"^+^ 



to construct {u"^~^^ , p"^^^] 



the sequence {{u'^,p^,'r]"^)}m=iJ which are recorded in the following theorem. 

Theorem 3.12. Assume that the initial data are given as in Section 3.3.1 and satisfy the 
{2NY^ compatibility conditions (3.83). Suppose that £,o < £o/2, where Eq is the smaller one 
from Theorem 3. 5 and Proposition 3. 8. There exist a < Mq < oo and a < Tq < 1 depending 
on So and k such that ifT<TQ, then for each m > 1 we have the estimates 



^{u"',p"^) + ^{rj"^) < Mo, and £.{7^) < Eq. 



(3.104) 



Proof. We claim that there exist a (large) < Mq < 00 and a (small) < Tq < 1 so that if 
T <To, then the following hold: if (u'",r/™) for m > are known and satisfy the estimates 

^2n{u"') + ^2Jv(r?"*) < Mo, (3.105) 

then there exists a unique triple ('u"*+^,p'"+-'^, 77"*+^) that solves (3.103) and obeys the estimates 

A{u'^+\p'^+'^) + .^(r/"*+^) < Mo, and £(17"^+^) < £o- (3.106) 

We remark that in the following determination Mq depends on Sq and To depends on Mq; so 
these two constants are determined by the initial data (and k). 

First, if we take Mq > P2{£q), then (3.105) holds for m = 0. Now assume that (3.105) holds 
for m > 0, we prove that (3.106) holds for m + 1 for appropriately chosen Mq and Tq. We 
first use {u^, rf^) to construct the 77"*+^ via Theorem 3.6. Moreover, by the estimate (3.55) of 
Theorem 3.6 together with the estimate (3.95) of Lemma 3.10 and (3.101), we have 



< P^{£q) + P3(l + d2N{vn)^2N{un- 



(3.107) 



We need an appropriate estimate of 52Ar(?7™') and ^2n{u"^), and the key point is to get a power 
of T in the following estimate. First, Lemma 2.6 allows us to bound, by (3.101) and (3.105), 



2N-1 



4N-2j 
2N 



2N 

+ 2j2\\div" 

j=0 



j^2ff4N-2j+l 



l^cx>ff4N-2j+l 



<P,{£o) + 2TY,\\div'' 

j=0 

< Pi(^o) + 2T^2n{v"') < PiiSo) + 2TMo. 



(3.108) 



On the other hand, using the usual Sobolev interpolation together with Lemma 2.6, (3.101) and 
(3.105) again, we obtain 



^2iv(u" 



2N-1 

3=0 



2N 



4iV-2j-l/4 



u 



j=0 



2N 





1/2 











£2jj4,N-2j+3/4, 

3/2 

J2JJ4N-23+1 



<P^{Sq) + 2Y,T^'^ 

j=0 

< Pi{£q) + CT'I^Si2N{un < Pi{£q) + CT^/^Mq. 

If we take T < min{^o/Mo, {Sq/Mq)^}, then we deduce from (3.108)-(3.109) that 

^2iv(r?"^) + ^2n{u'^) < Pi{£o) + £o< Pa{£q). 
Plugging the estimate (3.110) into (3.107), we obtain 

^(^m+i) < p^(^^) + p^^i + Pa{£o))Pa{£o) < P^iSo). 



(3.109) 



(3.110) 
(3.111) 
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This proves the bound of ^(77"*+^). Now, by again Lemma 2.6 and (3.111), we have 

- 2 II ' llL°°H4iV -t ll'-'t'/ ||^oo^4JV-l (3.112) 

- 2 ^ ' ' - 2 2 

If we restrict that T <£q/P^{£q), then we have 

£(r/™+^) < eo. (3.113) 

Now the bounds (3.111) and (3.105) and the smallness (3.113) allow us to use (ti™, ry™'"'"^) 
to construct the {vJ^'^^ ,p'^'^^) via Theorem 3.5. Moreover, by the estimates (3.50) of Theorem 
3.5, (3.101), the estimate (3.94) of Lemma 3.10, (3.111) and (3.105), we have 

^{u'^+\p^+^) < Pe{l + Mv""^^) + ^iv"'^^)) exp {Pe{l + div^'+^W) 

X (ihoIlL + UF\F') + ^{F\u^, r?-+i), F^{r,^+'))) 

< Pe{l + Pi{£o) + P5{£o)) exp (^6(1 + ^^5(^0))^) 

X (Pi(^o) + P7(l + ;?(r/™+')) i^ir,^^') + mN{u^)f)) ^^'^^^^ 

< Ps{l + £0) exp (P8(l + £o)T) (Pi(^o) + ^7(1 + P5{£o)) [P5{£o) + (^'4(^0))')) 
<P9(^o)exp (p8(1 + ^o)t) . 

If we restrict that T < l/Psil + £0), then 

A{u^+\p'^+^) < Pw{£o). (3.115) 
W can now fix Mq := P5{£o) + Pio(^o) and then take 



To := min < 



£0 1 £0 



m£o) + Plo{£o) ' V P5{£o) + ^'io(^o) / ' P5i£o) ' ^^8(1 + £0 



(3.116) 



Hence if < T < Tq, then we deduce the estimates (3.106) from (3.111), (3.113) and (3.115). 
We thus conclude the claim and complete the proof of the theorem by the induction. □ 

Now we can produce a unique local strong solution to the nonlinear K-problem (3.1) with 
K> 0, which is the main result of this section. 

Theorem 3.13. Let k>0. Assume that uq G H'^^{n),r)o G if^^+i(S) with £0 < eo/2 for the 
universal eo > from Theorem 3.12, and that the initial data dfu^O), etc are as constructed in 
Section 3.3.1 and satisfy the (2NY^ compatibility conditions (3.83). There exists a < Tq < 1 
such that there exists a unique solution triple {u,p,r]) to the problem (3.1) on the time interval 
[0, Tq] that achieves the initial data and satisfies 

A{u,p) + R{r]) < Pk (\\uo\\Ij^ + ||??o||4iv+i) O'lT'd ■C(^) < £0. (3.117) 

Proof. The uniform estimates oi ^{u"^^^ ,p"^'^^)+^{r]'^~^^) in (3.104) allows us to take weak and 
weak-* limits, up to extraction of a subsequence, so that the sequence (-u™", r/™") converges to 
a limit (u,p, rf) in the weak sense of the norms in the definition of ^{u,p) + ^{rj). If we recover 
the dependence of Mq and Tq on the initial data and n as Mq = Pk{£o) = -fK(l|^^o||4Ar+||??o||4Ar+i) 
and To = Tq, then we have from (3.104) that for < T < Tq, according to the weak and weak-* 
lower semicontinuity of the norms in A(u,p) + M.{r]), 

^{u,p)+^{r]) < P« (||uo||4^ + ||??o||4iv+i) • (3-118) 
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Similarly, the uniform smallness of £,{rj'^'^^) in (3.104) implies 

&{ri) < £0- (3.119) 

The convergence of a subsequence of {u"\p"\ ry™) does not guarantee that we can simply pass 
to the limit in the system (3.103) in order to produce the desired solution to (3.1). However, 
by restrict Tq to be smaller if necessary, the uniform estimates (3.104) will allow us to show 
that (u™' , p"^ , r]™') is contractive in certain norms of lower order regularity; the proof is very 
similar to Theorem 6.2 of [12] and we omit the details. This implies that the whole sequence 
of (ii™',p™', Ty™') converge strongly to the limit {u,p,T]). This strong convergence together with 
the higher order regularity weak convergence above allow us to interpolate to derive a strong 
convergence result that is more than sufficient for us to pass to the limit in the system (3.103) 
for each t G [0,ro^]. Then we find that the limit {u,p,r]) is a strong solution of the K-problem 
(3.1) on the time interval [0, TjJ'] that achieves the initial data. The uniqueness of the solution 
to (3.1) with K > satisfying (3.118) follows also from the contractive argument. □ 



4. Local in time theory 

The aim of this section is that first for each fixed a > we will pass to the limit as k — )• 
in the regularized /c-problem (3.1) to produce a local unique strong solution to the original 
problem (1.11), and then we will pass to the limit as cr —>■ in (l-H) to rigorously justify the 
zero surface tension limit within a local time interval. 



4.1. Preliminaries. To estimate (^-u, ^p) for j > 1, we will recall that {Dlu,dlp) solves 



f dt{Dlu) - Aa{DIu) + VAidip) = F^^i in ^ 

diY A{Diu) = inn 

SA{dip, Diu)N = F^'^ on S 

Diu = on Sb 



(4.1) 



in the strong sense for j = 1, . . . , 2N — 1, and in the weak sense of (3.5) for j = 2N. Here the 
vectors F^'^ and F'^'^ are as defined by (3.19) in which F^ and F'^ are as defined by (3.60). 

We will restate the estimates of F^'^ and F^'^ with a slight modification from Lemma 3.2. 
To this end, we define 



2N-1 

UF\F'):= \\DiF' 



]^2ff4N-2j-l 



+ 



2N 



+ Tmf 



(4.2) 



and 



^{u,p) ■= 



2 



2N-2 



2N-\ 



+ 



2 



i=0 
2N-2 



l,2ff4N-2e- 



u 



£=0 



l^ooJj4N-2l-2 



(4.3) 



J^2fj4N-2e-2 



+ 



dip 



l^oa }J4N -2t-3 



Remark 4.1. We can deduce an estimate of^{u,p) in terms ofT. Indeed, we apply Lemma 
2.6 to obtain that for i = 0, . . . ,2N - 1 



dfu 



< 



l^oo ]^4N -21-2 

and for £ = 0,...,2N -2 

2 



IP 



lj00fj4N-2l-3 



< 



afn(O) 



afp(o) 



4:N-2i-2 



+ 



dfu 



l2 ^4N-2l-2 



+ 



4Ar-2^-3 



+ 



%p 



l^2ff4N-2e-3 



+ 



dt'u 



l2fj4N-2i-2 



l2fj4N-2i-3 



(4.4) 



(4.5) 
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Then we have 



2N-1 



U 



2N-1 



l2[jiN-2i 



+ 



dip 



<€o{u,p)+T Wt 



u 



+ 



l2jj4N-2e-l 
2 



(4.6) 



dip 



^ < €o{u,p)+T^{u,p), 



l^ooJj4N-2e 

where ^o{u,p) and ^{u,p) are defined by (3.62) and (3.47), respectively. 

Wc now record a version of estimates of F^'^ and F^'^ , recalling 5(r/) as defined in (3.22), 



Lemma 4.2. It holds that 

2N-1 



-1 + \\F^'^ 



L2H4N-2j- 



(4.7) 



Proof. The estimate (4.7) follows by taking m = 2A^ - 1 in (3.24) and (3.26) of Lemma 3.2; the 
term ^(F^,F^) therein can be replaced by ^*{F^,F^) is due to the definition (3.19). □ 

We also record here a version of estimates for the difference between dtv and Dtv. 

Lemma 4.3. For j = 1, . . . ,2N - 1, 



dfu — Dfu 



- Di 



u 



d^u - dtD 



]^2fj4N-2j+l 
2 

2 



<P(l+;?(r/));?(n,0), 



2N-1 
t 



U 



L'>°H'> 



<p{i+d{vmiu,o), 



and 



+ 



<P{l + d{ri)) (d{u,0) + 



,2N-1„ 



u 



2 



|2 



(4.8) 
(4.9) 
(4.10) 

(4.11) 

□ 

(4.12) 



Proof. These estimates follow similarly as Lemma 4.2. 

Next, to estimate djrj for j > 1, wc will recall that 

dtf] = kAt] + + F^ on E. 

Finally, to estimate {u,p,r]) itself without temporal derivatives, we shall write the system 
(3.1) in the perturbed form: 

' dtu - Au + Vp = G^ inn 

div u = G"^ in Jl 

{pI-B){u))ez = {r]-aA^ri)e^ + G^ on S (4.13) 

dtrj — kA^t] — Us = + G^ on S 

u = Q on Sfc. 

Here we have written the vector G^ for 

G] = dtfjbKdsUi - UjAjkdkUi + {AjkAji - Sjkdji)dkiUi + AjkdkAjidiUi - {Aij - Sij)djp, (4.14) 
the function G^ for 

G'' = -{Aij - Sij)djUi, (4.15) 
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'{K-l)d3Ui+AKd3U3\ 

{K - l)d3U2 + BKd3U3 

2{K - 1)53«3 / 



(4.16) 



the vector for 

(p — T] — 2{diui — AKdsui) \ 
-d2Ui - diU2 + BKd^ui + AKd3U2 
-diU3 - KdsUi + AKd3U3 J 

(—d2U\ — d\U2 + BKd^ui + AKd3U2\ 
p - ri - 2{d2U2 - BKd3U2) 
-d2U3 - Kd3U2 + BKd^us J 

- adiv*(((l + |V*r/|2)-V2 _ l)V,r7)AA, 
and the function for 

G^ = -n-V*7?. (4.17) 

We want to apply the horizontal spatial derivatives to (4.13) with a G so that \a\ < AN. 
However, the lower boundary S;, may not be flat, so we are not free to apply such derivatives 
to (4.13). The idea then is to localize away from the lower boundary. To this end, we introduce 
the cutoff function x £ C'^(IK) so that 

suppx C [-36_/4, 1] and xixs) = 1 for X3 G [-b-/2, 1/2]. (4.18) 

Multiplying the equations in (4.13) by X; we find that {x'^i-XP^rj) solve 

'dt{xy) + V(XP) - A(xn) = xG^ + 

div(xti) = xG^ + H"^ mVL 

{{XP)I - Hxu))e-i = (r/ - c7A,j?)e3 + G^ on S (4.19) 

dtT] — nAitT] — [xus) = K'^ + on S 

X^i = on S(,, 



where 



= 53x(pe3 - 2d3u) - dlxu and H"^ = d3XU3. 



(4.20) 



Since now xu and XP have the support away from S;,, we could apply d"" (a G N^) to (4.19). 

4.2. Local well-posedness of (1.11) with surface tension. In this subsection we fix a > 0, 
and our aim is to pass to the limit as k — >■ in the K-problem (3.1). We define a special quantity 
for this subsection as 

/•^ , 2 

^'^ = A'^{u,p,r}) : = sup i£{t) 



0<t<T 



where 



and 



27V 



j=0 



2 

4JV-2j 



2N~l 



27V 



-2j ^ II * 4Ar-2j-l II ^ 



j=0 



27V 



j=0 



27V-1 



4JV-2j+l 



j=0 



j=Q 

2 

47V-2j 

27V+1 



4JV-2j+l 



(4.21) 



(4.22) 



(4.23) 



4JV-2j+2 



+ ll^ll47V+2 + ll'7ll47V+3/2 + E 

J = l 

For each k > 0, we let {u,p,r]) (dropping the K-dependence temporally) be the solution to 
the problem (3.1) on [OjTo^] produced by Theorem 3.13. Moreover, Theorem 3.13 implies that 
for < r < Tq^, we have .^'^ < 00 and the following smallness of rj: 



S,{ri) := sup_||77(i)||4jv_i/2 ^ 



0<t<T 



(4.24) 



for the universal constant £0 ^ 1- Our goal is to derive the ^-independent estimates of the 
solution on [0, Tq], and the key step is to derive the uniform bound of We will see that the 
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smallness (4.24) is crucial for closing the ^-independent energy estimates since we can not get 
a power of T in front of the estimates of some terms stemming from the mean curvature. 
We recall from the construction of the initial data in Section 3.3.1 that 



2Ar+i 



4Af-2j+3/2 



<f2V(0)<P(^o), (4.25) 



where £q := ||'Uo|l47v + II^IIItv+i' ^^"^ ^2N^^) defined by (1.12) with that o" = 1, n = 2N 
and t = 0. We also recall from the construction of ^ in Section 3.3.2 that, by (4.25), 



27V 
3=0 



2N 



L2H4N-2j 



2N 
j=0 



4N-2j+l 



< P{So) 



(4.26) 



4.2.1. The K-independent energy estimates with temporal derivatives. We first record an esti- 
mate of the nonlinear terms ^*{F^,F^) (defined by (4.2)) and an estimate of ^{F^) (as defined 
in (3.52)) in terms of A'' (defined by (4.21)). 

Lemma 4.4. Let F^{u,r]), F^{r)), F'^{u,ri) be defined by (3.60). Then 

UFHu,r,),F^{r,)) < P{S^) + P{T^'^^'') (4.27) 

and 

d{F\u, r,)) < P{£o) + P(rV4j^-). (4.28) 

Proof. We first prove (4.28). We recall from the estimate (3.95) of Lemma 3.10 that 

d{F\u, r?)) < P(l + d2N{v))d2N{u), (4.29) 

where ^2Niu), d2N{v) ^^^^ defined by (3.92) and (3.93), respectively. Wc then estimate d2Niv) 
and d2N{u). We recall from the estimates (3.108)-(3.109) that, by the estimates (4.25) and the 
definition oi ^'^ (4.21), 



and 



2N~1 

^2iv(r?) < E W^tm 

2N-1 

;?2Ar(u) < E 11^^(0) 

i=o 



2N 



4N-2j 



+ 2TY,Mv 



3=0 



]^caJj4N-2j+ 



CTV4^2Ar(u) < P{So) + CT^/^^". 



4Af-2j-l/4 

Plugging the estimates (4.30)-(4.31) into (4.29), then we have 

^{F\u,ri)) < P(l + PiSo) + 2W)(P(^o) + CtV4^«) < p^^^^ ^ P(tV4^«). 



(4.30) 



(4.31) 



(4.32) 



This gives (4.28). 

Now we turn to prove (4.27). We refine the estimate (3.94) from the proof of Lemma 3.10 as 



(4.33) 



Here the reason that wc are able to use 5^(r/) (as defined in (3.22)) to replace ^{r]) (defined 
by (3.53)) in (3.94) is due to that there is at least one temporal derivative in the definition of 
5^*(P^,P^); the appearance of ^{rj) in (3.94) is due to the estimate (3.99). We then appeal to 
an estimate of div)- First we trivially have 



II^IIl2//4]V+1/2 < r ||?7||^ooJ|/4]V+l/2 ^ TK'^. 

Next for j = 1, . . . , 2N, we use the Sobolev interpolation to obtain 

2 



]^2jj4N-2j+3/2 



< ^1/2 



Laofj4N-2j+l 



]^2fj4N-2j+2 



(4.34) 
(4.35) 
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Now for j = 2N + 1, we must resort to that rj is the solution to the regularized surface K-problem 
(3.51); then we have, by (4.35) and (4.28), 



(4.36) 



Hence, we conclude from (4.34)-(4.36), by further applying Lemma 2.6, 

2N ^ 2N+1 



m<\m\\lN+T.\rtm 



]^2ff4N-2j+3/2 



4N-2j+l/2 

< P{8o) + W + T^/'^a'' + PiT^I^^") < P{So) + PiT^'^A""). 
Plugging the estimates (4.31) and (4.37) into (4.33), we deduce 

d.{F\u,'n),F\r])) < P{l + P{£o)+P{T^'^ii'')) (p{£o) + PiT^'*^"") + {2T^'^ii''f) 
<P{So) + P{T^I^R''). 

This gives (4.27), and we then conclude the lemma. 

We shall now use the equations (4.1) to estimate (dfUjdlp) for j > 1. 
Proposition 4.5. Suppose that {u,p,r]) is the solution to (3.1). Then 



(4.37) 



(4.38) 



□ 



2N 



2N-1 



U 



L2H4N-2j+l 



+ 



g2N+l 



U 



{XtY 



L2H4N-2j 



2N 



2N-1 



l^ooJj4N-2j 



l^ooff4N-2j-l 



< (p(^o) + P(rV4j^«)) exp ((p(4) + P(ri/4j^'')) r) . 



(4.39) 



Proof. First, taking j = 2N — 1 in (4.1), we then apply the estimate (3.12) of Theorem 3.1 to 
find, using the estimate (4.7) of Lemma 4.2, 



+ 



2 



+ 



(Xt)* 



+ 



dtD, 

+ 



2N-1 



u 



L°°H0 



+ 



u 



2 



■'t 



P 



< P{1 + /C(r?)) exp (P(l + /C(??))T) (p(^o) + HP^'^^-^ 

, II p3,2Af-l||2 I 11^1 /pl.SAf-l p3,2Ar-U||2 \ 

+ \\F llL2/f3/2 + ||Ct(^ -F' )\\(Xt)*) 

P{l+^{v))exp{P{l+^{ri))T) (^P{So)+^.{F\F^) + ^{u,p)) := Z. 



(4.40) 



< 



We will allow Z to change in the sense that the polynomial P changes from line to line. 
While for j = 1, . . . , 2N — 2, we rewrite (4.1) as 

' -Aa{d{u) + Va{4p) = -dt{Diu) + F^'^' in Q 



div^(L»^u) = 
5^(a/p, Diu)J\f = F^'^ 
Diu = 



in i7 
on S 
on 



(4.41) 
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We apply the ^-Stokes elliptic regularity theory of Proposition 2.17 with k = 4iV and r = 
4N - 2j + 1 < AN - 1 to (4.41) to obtain, using the estimate (4.7) of Lemma 4.2, 



< 



< 



< 



2 



+ 



dip 



dt{Dlu) 



]^2H4N-2j-l 



L2H4N-2j 

^11 llL2H4^-2i-l ^ II llL2_ff4JV-2j-l/2 



l2fjiN-23-l 
2 



+ 



RDlu 



(4.42) 



L2HiN-2j- 



A simple induction on (4.42) yields 

2N-2 



< A 

W t 



2N~1 



l2ff4N-2j+l 
2 



+ 



dip 



u 



L2H'iN-2j 

1 Tp3\ 



(4.43) 



+ MF ^F-") + P(l + d{ri))^{u,p). 



Hence, we combine (4.40) and (4.43) to have 

2N-1 



u 



]^2ff4N-2j+l 
2 



+ 



dip 



+ 



dtD, 



2N-1 



l2H4N-2j 
2 



+ 



dfD 



2Af-l„ 



+ 



u 



2 



(4.44) 



u 



+ 



dr-'p 



< Z. 



Now wc want to derive more estimates based on the estimate (4.44). By the estimates (4.8)- 
(4.10) of Lemma 4.3 and (4.44), wc have 

2Af-l „ 2JV-1 „ 2JV-1 



l2jj4N-2j+l 



u 



]^2H4N-2j+l 



j=i " j=i 

<Z + P{l + ^{r,mu,0)<Z 



L2H4N-2j+l 



(4.45) 



and 



o2Ar-i 
9t u 



|2W„,||2 



t "IlLooffa < 



+ 



L°°H'^ 
2N-l„ ^ 



+ 



^2iV-l_^2Ar-l^ 



dtDf-'u 



+ 



2N-1 

t ' 



2 



(4.46) 



<Z + P(l+J(r/))i?(n,0) <Z. 
By the estimate (4.11) of Lemma 4.3, (4.44) and (4.46), we have 



\ t "■||l2^1 



Q2N+1 



< 



+ 



[XtY 



2N-1„ 
2 



<Z + P{l+d{r]))id{u,0)+ dt 



32Af-l 



2 



+ 



+ 



u 



2 



2 

{XtT 

< z. 



(4.47) 



On the other hand, we use Lemma 2.6 and (4.44)-(4.45) to deduce 

2Ar-2 



dip 



y + 

^-^ II J^^ll4N-2j 
2iV-2 2 

< y b/n(o) + 

II ^ ^ ' 4N-2j 

<p{io) + z<z. 



lj00ff4N-2j-l 
2 



2N-1 



dlp{0) 



4N-2j- 



-1+ E i^'^ 



j^2ff4N-2j+l 



+ 



dip 



(4.48) 



L2ff4N-2j 
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We can now conclude the estimates for {dfujOfp) for j > 1 from (4.44)-(4.48) as 



2N 



2N~1 



U 



L2H4N-2j+l 



+ 



{XtY 



L2H4N-2j 



2N 



u 



ljOoH4N-2j 



2N-1 

+ E W^tP 

i=i 



(4.49) 



l^ooJj4,N-23-l 



< Z. 



Turning back to estimate we employ the estimates (4.37), (4.27), (4.6) and (4.25) to deduce, 
since ^iu^'p) < 

Z = P{1 + div)) exp (P(l + d{v))T) [P{So) + MF^F"") + 5(«,p)) 

< P(l + P{So) + P(rV4^«)) exp (P(l + P{So) + P(TV4^«))r) 

X (^P{io) + p(rV4^«) + r.^«) (4.50) 

< (p(^o) + P(rV4j^-)) exp (p(i + + P(rV4j^«))r) 

< (p(^o) + P(tV4^'^)) exp ((p(^o) + P(rV4^«)) t) . 

Plugging the estimate (4.50) into (4.49), we then deduce (4.39). □ 

We then use the equation (4.12) to estimate dji] for j > 1. 
Proposition 4.6. Suppose that {u,p,ri) is the solution to (3.1). Then 

2N+1 ^ 2N ^ 

^ < PiSo) + P(rV4^-) + \\ri\\l2H4N+2 . (4.51) 



j^2ff4N-2j+2 



+EF 



l^aoj^4N-2j- 



Proof. We recall from (3.58)-(3.59) that, by (4.25)-(4.26), 

2N+1 „ 2N „ 



£ocfj4N-2j+l 



' II l2jj4N-2j+2 II ' 

2N 

^ ||??|li2^^4JV+2 + E 11^^(0) 
3=0 

< M\l2H^N+2 + P{£o) + ^{FHu, 7j)). 

Plugging the estimate (4.28) into (4.52), we obtain (4.51). 



4N-2j+l 



+ 5(P^(«,?7)) 



(4.52) 



□ 



4.2.2. The n-independent energy estimates without tem,poral derivatives. To simplify notation, 
we define the energy and dissipation without temporal derivatives by, where we temporally keep 
the dependence of cr. 



^ II \\^ I II II I II II I II II 

^ — iPlUw "T llf'lUw-i "I" \\m\AN "T ^ ll'^IUw+i 



and 



V 



_l_ 2 II ||2 
+ ^ ll^ll4Af+3/2 



+ ll^ll4Af-l/2 ■ 

We also define the "horizontal" energy and dissipation with localization as 

^ = ||V*o^(xtx)||o + ||??||4;v + <^ Ihll4iv+i and V = ||V*o^B(xu)||g + Ka ||??||47v+2 : 

where x is as defined in (4.18). 

We first record an estimate of the nonlinear terms G'\ defined by (4.14)-(4.17). 

Lemma 4.7. Let Eq be the small universal constant in (4.24). Then we have 



(4.53) 
(4.54) 

(4.55) 



\4N- 



-L. 4-\\ n-^ 4-\\ 

-l"^ll^ ll4Ar"^ll"^ ll4Ar-l/2 11^ ll4Ar-l/2 



<£oVP + P(l + e). 



(4.56) 
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Proof. We note that all terms are at least quadratic. Then we apply the spatial derivatives d" 
(for the appropriate a depending on G*) and expand using the Leibniz rule; each term in the 
resulting sum is also at least quadratic and can be written in the form XY or ZZ. Here Y 
is the factor involving the highest derivatives and X is a product of factors involving very low 

derivatives; Z is the factor involving fewer derivatives than Y so that the various norms of Z 
can be bounded by P(l + Wc then focus on how to estimate the terms of form XY. 

II 1 1|2 

We first estimate the ||q for a G N with \a\ < AN — 1. Note that those factors 

appearing as Y of d'^G^ are d'^V^u,d'^Vp,d°'dtfj,d'^V'^fj. Then we have 

II^^V^ullo < VV, ||a"Vp||o < VV, (4.57) 

ll^"^t^ll(oHl(f^))* ^ ll^*^ll4JV-l ^ ll^t^ll4iV-3/2 ^ (4-58) 

ll^"V'^ll(o/^Mn))* ^ ll^'^^ILiV-i ^ ll^ll4Ar+i/2 ^ (4-59) 



We can bound certain norms of X related to the last two factors by P{1 + ^)v However, to 
bound the X related to the first two factors we need to resort to the factor that such X only 
involves the spatial derivatives of fj (up to second order) ; hence we can bound the certain norms 
of such X by P(l + £{'rj))^/£^ < eq, using (4.24). In light of the analysis above and using 
that P(l + £)€ < P(l + ^), we conclude the estimate of as in Lemma 3.2. 

We now estimate the ||9"G^||^y2 a G with |a| < AN — 1. Note that the factors 
appearing as Y of d'^G^ are d'^p, 9"Vn, d°''Vfi and ad^V^r]. For the first three factors, we have 

l|5>ll/fi/2(E) ^ IblUiv < ll^"Vn||^i/2(s) < II^^IUat+i < VP- (4.60) 

||5°Vr?||^_i/2(s) < ||Vr7||^4^-i/2(s) < MUn+i ^ M\4n+i/2 ^ (4-61) 
To estimate the last factor we need to use the fact that there is a a. Hence we can bound by 

^ ||5"V^r/||_^/2 ^ ^ M\4N+3/2 < (4-62) 

Then arguing as for the G^ term, we conclude the estimate of G^. 

Finally, the estimates of G^ and G^ follows similarly as that for G^ and G^ and we then 
conclude the estimate (4.56). □ 

We now state the estimates of {u,p,'r]) itself without temporal derivatives as follows. 

Proposition 4.8. Suppose that {u,p,rj) is the solution to (3.1). Then 

S + [ V< P{£o) + P{1 + J^")Ti/2 + \\dtu\\l2H4N-, + \\dtp\\l2H4N-2 . (4.63) 
Jo 

Proof. We first derive the evolution of £. We claim that 



£+ f V < P{£o) + P{1 + ^'')T^/^ + eo [ V. 
Jo Jo 



(4.64) 



To prove the claim, applying 9° with q G so that |a| < AN to the first equation of (4.19) 
and then taking the dot product with (}"n, using the other equations, we obtain 

1 d_ 
2di 



[ r(xu)l'+ / \d"r,\' + a\d'^VM')+l I 
Jo. Js / ^ Jn 

+ / K |V*5"r?| V «a |A,5"r?p = / xd'^u ■ {xd'^G^ + d'^H^) (4.65) 
7s Jn 

+ [ xdy{xd''G^ + d''H^)+ f -df'u-&'G'^ + {&'r]-al^^d'^r]){Kdf'^ + &'G'^). 
Jn -/s 

We will estimate the terms on the right side of (4.65), beginning with the terms involving 
and H"^. Since x is only a function of X3, we have that 

d^H^ = dsxid^pes - Id'^d^u) - dlx&^u and d^H'^ = d^xd'^u^. (4.66) 



'3 11 

l4Ar-l/2 

(4.69) 
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This allows us to estimate 

/ x5"w • d^'H^ + xd'^pd'^H'^ < \\d^u\L {\\dy\L + \\d'^u\\^) + ||5>|L \\d'^u\L 
Jn (4.67) 

^ \\u\\iN {WpUn + II^^IUtv+i) < ^Vp. 

We now turn to estimate the terms involving G'\l < i < 4. Using Lemma 4.7 and the fact 
that 5" only involves the horizontal spatial derivatives, we may bound 

Jn 

< Wd'^uW, + II^VIIo \\d''G\ < Wd'^uW, \\G%^_^ + \\dy\\, \\G\^ (4.68) 

< VV (^eoVv + P(l + e)) < P(l + (E)VV + CeoV. 

Similarly, we use Lemma 4.7 along with the trace embedding to find that 

< (eoVv + P(l + €)) < P(l + €)VV + CeoV 

Now we estimate the term. For the first term for G^, we use Lemma 4.7 to have 

(4.70) 

< V€ {eoVv + P(l + €)) < P(l + €) + eoV. 
For the other term, we need to do the splitting: 

df^G^= Ga,pdN^r}-df^-^u + df^V^r}-u:=I + d''V^r]-u. (4.71) 

0</3<a 

Note that the highest derivatives appearing as F of / are d'^u and d°'~^V*ri for some ^ G 
with 1/3 1 = 1, and we have 

ll^"«ll/fV2(s) ^ ll^lUiV+i < 5"-^V*r? < ||r?||4^+i/2 < ^. (4.72) 

Then arguing as in Lemma 4.7, we may bound 

ll^ll 1/2 ^ ^ (so VP + P(l + €)) . (4.73) 

Hence, we have 

-aA^d'^ijI < a ||A*a°r/||_i/2 ll-f|li/2 ^ <^ ll^ll4JV+3/2 ll^lli/2 

< VV (^eoVv + P(l + (£)) < P(l + e)Vv + eqV. 
On the other hand, we integrate by parts to obtain, by Lemma 2.2, 

/ -aA^d''rid''V,v ■u = [ adid''r){V*did''r) ■ u + V*^"?? • diu) 
7s is 

= [ -a div* u ^^'^y^ + adid'^-qV^d^T] ■ diU (^•'^^) 
Js 2 

< ||V*m||j|^2(2) \W *d"-r]\\_ii2 *d^v\\i/2 ^ SV^. 
Hence, in light of (4.70), (4.74) and (4.75), we may bound 

j {d'^ri - aA^d'^ri)d"G'^ < P(l + (B)Vv + CeqV. (4.76) 



(4.74) 
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Finally, using Cauchy's inequality, we have 

/ (9"77 - aA,a"r?)K5"* < k ||a"r/||^ + ]-Ka ||A,a"r/||^ + C ||*||^^ . (4.77) 
Now, in light of (4.67)-(4.69), (4.76) and (4.77), we deduce from (4.65) that 



(4.78) 



|2 



(4.79) 



< P(l + £)\/p + Ceo^ + C 11*1 
for all a € with |a| < AN. We then integrate (4.78) in time to have, by (4.26), 

£+ fv<P{£o)+ f P{l + '^){\ + ^/V) + e^{D + \\nlN 
Jo Jo 

< P{So) + P{1 + A^)T^/^ + so f V. 

Jo 

This completes the proof of the claim (4.64). 

Now we will show that T) is comparable to T> so that we can improve the energy evolution 
obtained in (4.64). We will employ the elliptic theory to improve the estimates of {u,p). Notice 
that involves the high spatial derivative of r/, say, V^r], neither £ nor V contain an enough 
K-independent regularity estimate of r] that allows us to apply Lemma 2.15. But it is crucial to 
observe that we can get higher regularity estimates of u on the boundary E from T>. Indeed, 
since S is flat, we may use the definition of Sobolev norm on S, the trace theory and Korn's 
inequality of Lemma 2.5 to obtain, by the definition of P, 

II l|2 _ II ||2 ^ II ||2 I ||v74Af/ \||2 

||lt||_ff4JV+l/2(s) — \\XU\\fj4N+l/2(^Y:) IIX^'IIl2(S) + 11^* VX'"'J||_H-1/2(S) 

<\\xu\\l + \\Vt''ixu)\\l<V. 

This motivates us to use Lemma 2.16. Observe that {u,p) solves the problem 

' -Au + Vp = -dtu + in 

div u = G'^ inU 

u = u on S 

u = on Hh- 



(4.80) 



(4.81) 



We then apply Lemma 2.16 with r = 4A^ + 1 to have, by Lemma 4.7, 



2 2 2 2 lllll^ Il2ll^ 2 

Mm+i + WpWan ~ Ibllo + ll^t«ll4Ar-i + \\G W^^-i + \\m + Ibllif4iv+i/2(s) 

|2 



(4.82) 



(4.83) 



< P(l + e) + eoV + \\dtu\\ij,_^ + V. 
Combining (4.64) and (4.82), we obtain 

£+1 \MIn^i + \\p\\In + ^^<^MIn^2^P{£o) + P{'^+^'')T^'^ + ^0 ! 'D. 
Jo Jo 

We now improve the estimates of we use the boundary condition 

r) - aA^r) =p- 2d^U2, - Ol on S. (4.84) 

We use the standard elliptic theory on (4.84) to obtain 

rt ft 

I 2||||2 ||||2 /III|2 1 1 2 1 1 2 1 1 2 

ll'?ll4Ar+3/2 + Nll4Af-l/2 ^ ||p|| jj4iv-i/2(s) + II V^II//4JV-i/2(s) + ||G ||4^_;i/2 

< r lbll4Ar + lbll4Ar+i + eoV + P{1 + £) (4.85) 

^0 

<P(1 + .^«)tV2 + f \lp\ll^ + \\u\\Ij^^, + eoV. 
Jo 
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Now, in light of (4.83) and (4.85), we have 

S+ f V < P{£o) + Pil + Si^)T'^/^ + eo [ V. (4.86) 
JO Jo 

If wc have made £q be smaller (wc can achieve this by assuming £o be smaller in Theorem 3.13), 

then we can conclude (4.63) by further applying Lemma 2.6 to control the energy. □ 

4.2.3. The limit as n 0; local well-posedness of (1.11) with a > 0. Now we can produce a 
unique local strong solution to the original problem (1.11) for each fixed cr > 0, which is the 
main result of this subsection. In this subsubsection, we will allow the generic polynomial P to 
depend on a. 

Theorem 4.9. Let a > 0. Assume that uq £ H^^{n),rio G if'^^+^(S) with Qq < eo/2 for the 
universal eo > from Theorem 3.13, and that the initial data d^u{0), etc are as constructed in 
Section 3.3.1 and satisfy the {2NY^ compatibility conditions (3.83). There exists a < Tq < I 
such that there exists a unique solution triple {u,p,ri) to the problem (1.11) on the time interval 
[0, Tq'] that achieves the initial data and satisfies 

><^2Niu, P, v) < Pa (\\uo\\In + 11^011^+1) and&{rj)<£o, (4.87) 

where the energy functional lC\jq is the one setting a = 1 in (1.16). Moreover, 77 is such that 
the mapping $(-,t), defined by (1.8), is a (74iV-3/2 ^iff^QjyiQ'ppfiism for each t G [OjTq^]. 

Proof. A suitable linear combination of the estimate (4.39) of Proposition 4.5, the estimate 
(4.51) of Proposition 4.6 and the estimate (4.63) of Proposition 4.8 yields 

< (^'(^o) + P(ri/4.^«)) exp ((p(^o) + P(rV4^«)) r) 

+ P{So) + P(l + ^'')T^/^ + P(rV4^'^) (4.88) 

< (Pi(^o) + P(l + J^")Ti/^) exp ((p(fo) + p(rV4i^'')) r) . 

Note that here we have transferred the a dependence of £ and T> to the polynomials P, which 
we have allowed to depend on a. We name the polynomial Pi so that it can be referred later. 

We can now use a standard continuity argument (cf. Section 9 of [10] for instance) to infer 
from (4.88) that there exists a Tq" > that depends on the initial data and a but does not 
depend on k so that for < T < Tq', recovering the dependence of the solution on k, 

J^''(n^p^r/^) <2Pi(fo). (4.89) 

This K-independent estimates in turn imply that (u'^ , p*^ , rj'^) is indeed a solution of (3.1) on 
the time interval [0,rQ'] and yield a strong convergence of {u'^,p'^,rj'^) to a limit {u,p,r]), up to 
extraction of a subsequence, which is more than sufficient for us to pass to the limit as k — t- 
in the system (3.1) for each t € [0,Tq"]. We then find that {u,p,r]) is a strong solution of (1.11) 
on [0, Tq] and satisfies the estimates 

Si%u,p,r]) < P„ {\\u4l^ + ||%||4;v+i) (4-90) 

if we recover the dependence of Pi on o", where the quantity ^ is the one setting k = in 
(4.21). Moreover, if we restrict T^ to be smaller if necessary, then we have 

^{n) < £0 + \H\l2mN + WdtriWl^H^M-i < y + Ti^° < £0- (4.91) 

The uniqueness of the solution to (1.11) with o" > satisfying (4.90) follows as Theorem 3.13. 

To complete the theorem, we need to improve the estimates (4.90) to be (4.87). Comparing 
the definitions of ^ (as defined by (4.21)) and }C\^ (as defined by (1.16)), it means to improve 
the estimates of rj with temporal derivatives. For j = 1, . . . ,2iV + 1, we use the transport 
equation dtrj = F'^{u, rf) to have 



2 



^2jj4JV-2j+6/2 



dir] = di-^F^ = di'^F^ . (4.92) 



i2jf4Ar-2j+B/2 



2 



icx.^i-4JV-20-l)+l/2 
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But as in Lemma 3.2 we easily have 



]^2ff4N-2{j-l)+l/2 



These two imply, by (4.90), 

2N+1 



E 



\div 



]^2ff4N-2j+5/2 



< ^0 ^ (^0)2 < j|||^^||2^ ^ ||r/o||4Ar+l) • 



(4.93) 



(4.94) 



This completes the improvement of the L^H^ estimates of dlrj for j = 1,...,2N + 1. We 
may further apply Lemma 2.6 to improve the L°°H^ estimates of &lrj for j = 1, . . . , 2N. For 

D'i^u ■ J\f. Then by Lemma 2.11 and 



J = 2N + 1, we shall use the geometric fact: d^^~^^r] 
since div^(£)|^n) = 0, we have 

2 



L°o^i--l/2 



lorn 



u\ 



2 ^ p f\\ ||2 1 II ||2 

^ ^0- ^^IPOIUat + ||%|l4iV+i 



(4.95) 

We thus conclude the estimate (4.87) in view of the analysis above. The C^^~^^^ diffeomorphism 
of $(-,t) for each t G [0, Tq'] follows from the regularity of rj and the smallness (4.91). □ 

4.3. Local zero surface tension limit of (1.11). In this subsection our aim is to pass to the 
limit as o" — > in the problem (1.11) within a local time interval, where cr > is the surface 
tension coefficient. We define a special quantity for this subsection as 

A" = A''{u,p,ri) := sup (£"(1)+ [ ^''{t)dt + 



0<t<T 



Jo 



where 



2N 



2N-1 



2N 



3=0 



and 



2N 



j=0 



. „ li * 4Ar-2j II * 4Ar-2.7-l 'll4JV+i 



2j- 



2N-1 



3=0 



4Ar-2j+l/2 



(4.96) 



(4.97) 



3=0 



u 



4N-2j+l 



j=o 



2 

4N-2j 
2N+1 



(4.98) 



4iV-2i+3/2 



+ ''■^ ll'?ll47V+3/2 + ll'?ll4Af+l/2 + E 

3=1 

We note that by the definition, 

^'^ = ^{u,p) +^{rj) + \\rj\\j^^ij4,N+i/2 +cr \\rj\\ 

+ ll^llL2^4Ar+3/2 , (4.99) 

where .^(n,p), 5^(??) arc defined by (3.47) and (3.22), respectively. 

For each a > 0, we let {u,p,r]) (dropping the a-dependence temporally) be the solution to 
(1.11) on [0, r^] produced by Theorem 4.9. Moreover, Theorem 4.9 implies that for < T < T^, 
we have .^'^ < oo and the following smallness of rj: 



S,{ri) := sup_||?7(f)||4^_i/2 < ^0 



0<t<T 



(4.100) 



for the universal constant <^ 1. Our goal is to derive the a-independent estimates of the 
solution on [0, Tq^], and the key step is to derive the uniform bound of A"^. The smallness (4.100) 
is crucial for closing the cr-independent energy estimates again since we can not get a power of 

T in front of some terms stemming from the mean curvature. 

We recall from the construction of the initial data in Section 3.3.1 that 

C:"(0) < /C^^(O) := £^^{0) + J'2n{0) < P{£^). (4.101) 

Here we recall that Smi^)^ ^2iv(0) and /C^jv(O) are as defined by (1.14), (1.12), (1.15) and 
(1.16), respectively. 



42 



ZHONG TAN AND YANJIN WANG 



4.3.1. The a -independent transport estimates. We shall first use the transport equation to derive 
the estimates for rj: 

dtT] + u ■ V*r/ = U3 in S 



Proposition 4.10. We have the following estimates for rj: 

h\\looH^N+„2 < exp (crV^v^) {\\vo\\In+i/2 + TA'') 

and 

5(r?)<P(a) + TV4p(^^). 
Proof. Basing on the transport theory of Lemma 2.19, we may obtain 

||??||^cx.f^4Jv+i/2 < exp lh(0lli?4Jv+i/2(s)C^i^ ^V^ + y Ih3(i)||j/4jv+i/2(s) c?*^ ■ (4.105) 

The trace theory and the Cauchy-Schwarz inequality implies 



r\Hm 

Jo 



i?4JV+l/2(s) 



dt< 



Tikwii 

Jo 



Then the estimate (4.103) easily follows from (4.105) and (4.106). 
To estimate S{r)), we first trivially have 

||'7|Il2^4JV+1/2 ^ ||r7||^c>o/f4iV+l/2 ^ T^. 

For _7 = 1, . . . , 2 AT + 1, we use the transport equation directly as dtr] = u ■ M to have 



12JJ4N-23+3/2 



l2fj4N-2j+3/2 



di-\u-M) 



1^2h4N-2U-1)-1/2 



By the trace theory and the Sobolev interpolation, we obtain that for £ = 0, . . . , 2A^, 



dim 



< 



L2HiN-2t-l/2^Y.) 



l,2ff4N-2e+3/4 



< 7^1/4 



1/2 

£j<x>ff4,N-2e 



3/2 

]^2ff4N-2e+l 



Trivially again, for £ = 0, ... , 2iV, 

2 



dlr] 



]^2jj4N-2t+l/2 



< T 



1^00 u AN -21+1/2 



Hence, we deduce from (4.109) and (4.110) that for j = 1, . . . , 2iV - 1 



0-1 



< 



^2j:^4iV-2(j-l)-l/2 



^^=0 



+ 



i2^4]V-2f-l/2(s) 

This together with (4.108) implies 



(4.102) 

(4.103) 
(4.104) 



(4.106) 

(4.107) 
(4.108) 

(4.109) 
(4.110) 



£^2j:^4JV-2£+l/2 



(4.111) 



(4.112) 



^2jj4JV-23+3/2 

We thus conclude (4.104) by further applying Lemma 2.6 to control the L^H'' norms in ^{rj). □ 

4.3.2. The remaining a-independent estimates. To derive the remaining cr-independent esti- 
mates for the solution, our strategy follows Sections 4.2.1-4.2.2. We first record an estimate of 
the nonlinear terms ^^^{F^ , F"^) (defined by (4.2)) in terms of A'^ (defined by (4.96)). 

Lemma 4.11. Let F^{u,r]),F^{r]) be defined by (3.60). Then 

UF\'^,v),FHv)) < P{S^) + T^I^P{R-). (4.113) 
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Proof. We also recall the estimate (4.33) from the proof of Lemma 4.4 that 

:S.{F\u,r^),F\r^)) < P{l+m) W + (^2iv(n))') (4.114) 
We recall from the estimate (4.31) from the proof of Lemma 4.4 that, since M2n{u) < 

2N-1 



Mu) < \\dlu{0) 

j=0 



(4.115) 



4Ar-2j-l/4 

Plugging the estimates (4.115) and (4.104) into (4.114), we have 

^.{F\u, ri),F\ri)) < P{1 + P{£^) + T^/^P{ir)) (p(^o") + rV4p(j^-) + + 2T^/''R''f 

<P{S^)+T^I^P{^''). 

(4.116) 

We thus conclude the lemma. □ 

For the nonlinear terms (defined by (4.14)-(4.17)), we have the following estimates. 
Lemma 4.12. Let Eq he the small universal constant in (4.24). Then we have 

\\gXn-i + \\G'\In + l|Gi4iV-i/2 + \\gXn-i/2 ^ eoVv + P{l + ^ 
Proof. The estimate follows in the same way as Lemma 4.7. 
Proposition 4.13. Suppose that {u,p,ri) is the solution to (1.11). Then 

^{U,p) + a ||ry||^ooj:^4JV+l + Cr^ ||??||^2^4JV+3/2 

< (^PiSS) + P(l + J^'")ri/4) exp ((^(^0"") + T^^^P{^")) ■ 

Proof. We recall the estimates for {dlujdfp) for j > 1 from (4.49) that, using the estimates 
(4.113) of Lemma 4.11, (4.101) and (4.6), 



(4.117) 
□ 



(4.118) 



27V 



2N-1 



U 



l2ff4N-2j+l 



+ 



2N 



[Xt) 

2N-1 



Looff4N-2j-l 



(4.119) 



< P(l + '^{n)) exp (P(l + 'S{ri))T) {P{£^) + MF\F^) + S{u,p)) 

< p{i + p{£S) + rV4p(^<^)) exp (p(i + P{£S) + rV4p(^^))r) 

< (^P{£S) + TV4p(^<^)j exp (p(i + £S + ri/4p(.^^))r) 

< (^P{SS) + rV4p(^<^)) exp ((p(fo°^) + rV4p(^<^)) r) . 

In view of the estimate (4.117) of Lemma 4.12 and Proposition 4.8 (setting k = 0), we obtain 



+ 



\]^2jj4N + 11^11 ^ooj:^4iV + 



lLo°_f/4JV-l 



"I" ll^ll L2//4JV+3/2 + ||?/|| ^2£f4iV-l/2 + O' 1 1 1 1 j^oo Jif 4JV+1 + 1 1 ?7 1 1 j;^oo J|^4JV 
< PiSS) + P(l + J?-)Tl/2 + \\dtu\\l2H4N-^ + \\dtp\\l2H4N-2 . 

Suitably summing up the estimate (4.119)-(4.120) , we get 

^{U,p) + (T^ ||?7||^2jLf4JV+3/2 + ||??||^2£f4iV-l/2 + CT 1 1 1 1 /^oo £f4JV+l + 1 1 1 1 ^oo £f4iV 

< (^P{£S) + T^^^p{sr)'^ exp ((p(£:o") + r^/^p(j?'^)) t) + P{£S) + P(l + J^")Ti/2 
(^PisS) + P(i + .^^)rV4) exp (^(^Piss) + rV4p(^^)) r) . 



(4.120) 



(4.121) 



< 
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We thus conclude the proposition. □ 

4.3.3. The limit as a 0: the proof of Theorem 1.1. We can now prove a ^-independent local 
well-posedness of (1.11) which allow us to pass to the limit as a —> in the local time interval. 

Proof of Theorem 1.1. Assume that the initial data Uq and r]^ satisfy the assumptions of Theo- 
rem 1.1 for < o" < 1. For each cr, wc let {u"' ,p'^ ,r]'^) be the solution to (1.11) on [0, Tq"] produced 
by Theorem 4.9. By the estimate (4.103)-(4.104) of Proposition 4.10 and the estimate (4.118) 
of Proposition 4.13, in light of (4.99), we conclude that 

^ (Pi^S) + ^(1 + -t^")r^/^) exp ((p(^o") + rV4p(j^'^)) r) 

+ exp (ctV^VF) (ll%ll4iv+i/2 + T^'") (4-122) 

< [Piiss) + P(i + .^^)ri/4) exp ((p(£:o") + rV4p(^^)) r) . 

Here we name the polynomial Pi so that it can be referred later. 

Since (""o'^o) ~^ ^oi'^o^Vo) as cr — >■ 0, we can now use a standard continuity argument to 
infer from (4.122) that there exists a Tq > that depends only on the initial data but does not 
depend a so that for < T < Tq 

i^'^(«^p^r/-) < 2Pi(£:o-) = P{£^). (4.123) 

This cj-independent estimates in turn imply that (u"' ,p"', t]^) is indeed a solution of (1.11) on the 
time interval [0, Tq]. Moreover, we can improve the estimate (4.123) to be (1-17) by improving 
the estimates of rj with temporal derivatives as in Theorem 4.9. 

The (T-independent estimates (1.17) yield a strong convergence of (n'^,p'^, r/^) to a limit 
{u,p,r]) in any functional space that /C2^ can compactly embed into, up to extraction of a 
subsequence, which is more than sufficient for us to pass to the limit as cr ^ in (1.11) for each 
t G [0,ro]. We find that {u,p,r]) is a strong solution of (1.11) for cr = on [0,To] and satisfies 
the estimates (1.18). It can be shown as in Theorem 3.13 that the solution of (1.11) with a = 
satisfying (1.18) is unique. This in turn implies that indeed the whole family of the solutions 
{u'^ , p'^ , T)^ ) converges to the unique limit {u,p,r]). We thus conclude the zero surface tension 
limit within the local time interval [0,ro] and complete the proof of Theorem 1.1. □ 



5. Global in time theory 

The aim of this section is to prove a cr-independent global well-posedness of the problem (1.11) 
and hence justify the global-in-time zero surface tension limit in (1.11). Note that Theorem 1.1 
in particular implies that if £q < 5q <^ 1 {Sq will be determined later but is smaller than eo/2 
in Theorem 1.1), then there exists a unique solution {u,p,ri) (dropping the a dependence) to 
(1.11) on [0,To] for Tq = Tq{£q) > 0. Moreover, the solution satisfies the estimates 

^2^(7^0) < JC^2n{To) < PiSo) < So- (5.1) 

Hence loi 6 > (sufficiently small that will be determined later), if £q is sufficiently small, then 

Qmin) < S. (5.2) 

This smallness will be used to derive the global energy estimates throughout the rest of this 
section. Since tjq satisfies the zero average condition (1.19), we have 

J r]{t) = for i > 0. (5.3) 

Indeed, from the second and fourth equations of (1.11), we obtain 

j^J^V = J^dtr] = J^u-M = J^dwAU = 0. (5.4) 

The zero average (5.3) will allow us to use Poincare's inequality on S. 
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5.1. Preliminaries. To derive the global energy evolution of the pure temporal derivatives of 
the solution to (1-11) we shall use the following geometric formulation. Applying the temporal 
differential operator df for j = 0, . . . , 2N to (1.11), we find that 



dt{d{u) - dtfjbKdsidiu) + u ■ V A{diu) + div^ S^idip, 4^) = F^'^ in Q 

div^(a^n) = in Q 

S^idip, diu)J\f = {dir] - aA^dir])J\f + F^'^ on S 

dtidiri) = 4u-M + F^'^ on S 

dfu = on Sft, 



(5.5) 



where 



o<e<j 



(5.6) 



+dfAikdi ^dkiAimdmUi+AimdmUi)-df{uiAik)di ^dkUi-dfAikdi ^dkP^ , 

= - E C'j4Aikdl-%ui, (5.7) 
o<e<j 



o<e<j 



- a ^ C^jdfMdi-^A.r, - adi (div*(((l + |V*r?p)-V2 _ i)v,,,)Ar) , 

o<e<j 



and 



(5.8) 



(5.9) 



o<e<j 

We shall present the estimates of these nonlinear terms F^'^ , F"^'^ , F^'^ and F^'^ in the fol 
lowing lemma, at both 2N and N + 2 levels. 

Lemma 5.1. Let F^'^ be defined by (5.6)-(5.9), then the following estimates hold. 

(1) For 0<j< 2N, we have 

and 

(2) ForO<j<N + 2, we have 

ll^'llo + ll^'llo + ll^*(-^^''')|lo + ll^'llo + ll^'llo ^ ^2^+2 

and 



(5.10) 
(5.11) 

(5.12) 



llo ~ ^2N^N+2- 



(5.13) 



Proof. All these estimates, except the estimates of the o"-related terms in the definition (5.8) 
of F^'^ , are recorded in Theorems 4.1-4.2 of [13]. The highest derivative appearing in these 



cr-rclatcd terms is djVlr]. Since 



wc can deduce (5.10) and (5.12). □ 



To derive the global energy evolution of the mixed time-horizontal derivatives of the solution 
to (1.11) we shall use the following perturbed formulation, i.e., setting k = in (4.13), 

' dtu- Au + Vp = G^ in Q 

div u = G'^ in O 

{pI-I]){u))e3 = {r]-aA^r])e3 + G^ on S (5.14) 

dtT) — Us = G^ on S 

u = on Hi,. 
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We shall again need also the localized equations of (5.14), i.e., setting k = in (4.19), 

' dtixu) - A{xu) + V(xp) = xG^ + in 

div(xit) = + H'^ in 

((XP)/ - ID)(x^x))e3 = (r? - aA,r/)e3 + G3 on E (5.15) 

dtV - ixuz) = G^ on S 

Xii = on E;,, 

where we recall x ^^d H^,H^ from (4.18) and (4.20), respectively. 

We shall present the estimates of the nonlinear terms , , G^ and G^ in the following 
lemma, at both 2N and N + 2 levels. For this, we define one more specialized energy term by 



/C := 



Note that /C < £n+2- 

Lemma 5.2. Let be defined by (4.14) -(4.17). Then 
(1) It holds that 



(5.16) 



+ 



+ 



V 



4Af-2^3 



*0 



1/2 



+ 



V 



47V- 1^4 
*0 



(5.17) 



< 



1/2 



(^ZN")^ + ^J^2N, 



and 



V: 



^ 



+ 



V 



4Ar-1^2 




+ 



V 



4Af-l 
*0 



G' 



2 

1/2 



+ 



V 



4Ar-1^4 
*0 ^ 



1/2 



+ V 



r-%G%^ + a' 



IVf^G^I 



(5.18) 



1 1/2 



< ^2Af^2iV + ^^2N- 



(2) It holds that 



j2{N+2)-2^1 







+ 



72(Ar+2)-2^2 



+ 



^^2(iV+2)-2^3 



1/2 



+ 



^ 2(Af+2)-1^4 



(5.19) 



1/2 



< OCT OCT 



and 



72(Ar+2)-l^l 







+ 







+ 



V 



2(Ar+2)-1^2 







+ 



^^2(iV+2)-1^3 



1/2 



V 



2(Ar+2)^4 



*0 



(5.20) 



1/2 



'^2N^N+2- 



Proof. These estimates arc recorded with the slight differences in Theorems 3.1-3.2 of [13]. 
These differences result from the estimates of the G^ and G^ terms. For the G^ term, the 
difference results from estimating the a-related terms in the definition (4.16) of G^. The highest 
derivative appearing in these cr-related terms is crV^ry. We need to make use of this a factor 
so that these tj-related terms can be bounded in terms of £^ and P^. For example, to bound 

|2 

ll/2' 



|V47V-1(^3||^ we estimate a llVf^-^V^ 



^lll/2 - ^ ll^ll4Ar+3/2 



— v^2N- Hence, we can deduce 
the estimates of G'^ in (5.17)-(5.20). For the G^ term, the difference results from that we control 
the higher regularities of G^ than those in Theorems 3.1- 3.2 of [13], but the proof is similar. □ 



5.2. Global energy evolution. In this subsection, we will derive the global energy evolution 
of the solution to (1.11). We define the "horizontal" energies and dissipations with localization 
as follows. We define those only involving temporal derivatives by 



£^'' = ^\\VJdiu ^ + T.\\^tri , + 

j=o j=o j=o 



and Vi 



j=0 



u 



(5.21) 
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We define those locaHzed away from the lower boundary S;, involving both the horizontal spatial 
derivatives and temporal derivatives by 

n— 1 „ n—1 „ 



£^'^ = ||v,r-^(x«)||o + ||v*vi"-nx^)||o + E 11^' 

j=0 



2n-2j 



j=0 



2n-2j+l 



and 

We now define 



= ||Vi"-^B(xn)||^ + ||V.Vi«-iD(x«)||^ 
C = C'* + and P„ = + Vt. 



(5.22) 



(5.23) 
(5.24) 



Remark 5.3. Note that we only consider the energy evolution of localized terms away from E(,. 
The method employed in [13] involves another energy and dissipation pair for terms localized 
near E;,. Here, our modification of the method of [13] frees us from the need to introduce such 
a lower localization. 

We first derive the global energy evolution of the pure temporal derivatives of the solution, 
that is, the evolution of £n^. We first present the result at the 2N level. 



Proposition 5.4. It holds that 



d_ 
di 



and 



(5.25) 



(5.26) 



Proof. Taking the dot product of the first equation of (5.5) with dju and then integrating by 
parts, using the second, third and fifth equations, we obtain 

2 



2dtJn 



2 1/" 



= [ J{diu-F^'^ + dipdwA{diu))- [ Sik{dip,diu)Mkdiui 
Jn JT, 

= I J{diu ■ F^J + dipF^'^) - [ {dirj - aA^diriW ' 4^ + F^'^ ■ 4^- 
Jn Jt, 

The fourth equation further implies 
- / id{r,-aA,d{r,W-diu = - [ (Sf r/ - aA,a^r/)(5t(af r?) - F^-^) 



(5.27) 



f 

~~2Jth 



+ 



Jjdlri-aA,dir,)F''^. 



(5.28) 



Hence, we have 



+ 



+ 0- 



(5.29) 



= [ J{d{u ■ F^'^ + dipF"^'^) + [ -d{u ■ F^'^ + {d{ri - aA,diri)F^'^. 
Jn Jt. 

We now estimate the right hand side of (5.29). For the F^'^ term, by (5.10), we may bound 

|P''1„ £ V^V^A- (5-30) 



/ 

Jn 



Jdlu-F^'^ < ||J||^oo 



For the F^'^ and F^'^ terms, by (5.10) and the trace theory, we have 



J -dlu ■ F^^^ + {dlr] - (jA^&lri)F^^^ < 



d{i 



|_p4,j| 



< 



( 



+ 



(5.31) 
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For the F'^'^ term, we must consider the case j < 2N and j = 2N separately. For j < 2N, 
by (5.10) we have 

1^ JdfpF''^ < \\J\\^^ \\dip\l \\F''% < ^fvf^sj£',A- (5-32) 

The case j = 2N is much more involved since we can not control df^p. We are then forced to 
integrate by parts in time: 



[ Jd^pF^^-''' = [ 9t(J52^- V''^) - / d^'^-'pdtiJF 
Jn Jci Jci 

By (5.10), we may bound 

~ J^dr-'pdtiJF') < \\dr-'pl \\dt{jF')i < Vn^y^ 



Now we combine (5.30)-(5.34) to deduce from (5.29) that, summing over j, 

2N „ 



d 



dt 



a 



j=Q 



a 



2N^2N- 



We follow the estimates (4.17)-(4.21) of [13] to have 



J 



d{u 



> 



C\l £2n'^2N- 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



We then deduce (5.25) from (5.35)-(5.36). On the other hand, we may use (5.11) to obtain 



/ 

Jo. 



i2Ar-l„T72,2Ar 



jdf-'pF 



< \\J\ 



L°° 



This implies (5.26) and wc thus conclude the proposition. 

We then record a similar result at the N + 2 level. 
Proposition 5.5. It holds that 



(5.37) 

□ 



dt ^^^+2 



and 



- 2 [ jaf + V'"^- 



N+2 



(5.38) 



(5.39) 



Proof. The proof is similar to Proposition 5.4, except using (5.12)-(5.13) in place of (5.10)- 
(5.11). We thus omit the details. □ 

We now derive the global energy evolution of 82^^. 
Proposition 5.6. For any e G (0, 1) there exists a constant C{£) > such that 



jfii + '^tN ^ V^'^2N + V'^2n'^^2N + eV-2N + CisM^. 



(5.40) 



Proof. Applying 5" with a G N"^"'"^ so that oq < 2N — 1 and |a| < 4A^ to the first equation of 
(4.19) and then taking the dot product with d'^u, using the other equations as in Proposition 
5.4, we find that 



Jci Jn 



(5.41) 
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We will estimate the terms on the right side of (5.41), beginning with the terms involving 
and H'^. By the expression of and H"^ (4.20), we have 

/ Xd"u ■ dPH^ + xdyd'^H^ < ||5"u|L (||5>|L + \\d'^u\\^) + ||5>|L \\d'^u\L 
Jn (5.42) 

< \\d-u\\, iWdyi + \\d^u\\,) < K^u\\,^_,^^ 

We use the standard Sobolev interpolation to obtain 

||9"°^||4iV-2.o < l|9"°^llol|9"°^|iy_2.o+i ^ C^mf^'i-^m)^'-'^^', (5-43) 
where 6 = {AN — 2ao + 1)^^ € (0, 1). This and (5.42) together with Young's inequality imply 

Xd^ ■ d-H' + xdyd-H^ < mNf\l^2N)'-'^' < eV-2N + ^(e)^^. (5.44) 



/ 



We now turn to estimate the terms involving G', 1 < i < 4. Assume initially that \a\ < 4N—1. 
Then by (5.18) along with the trace theory, we have 



Jn 



(5.45) 



and 



(5.46) 



Now assume that \a.\ = 4N. Since ao < 2N — 1, we have that ai + a2 > 2, then we may write 
a = /3 + (a — /3) for some /3 G with \P\ = 1. Since |a — /3| = 4N — 1, we can then integrate 
by parts and use (5.18) to have 



X 





< 






Qa-HQl 














(5.47) 



< ||V 



For the G^ term wc do not need to integrate by parts: 



/: 



X^dyd'^G^ 



< WpWo ||v 



4Af-1^2 



For the G^ term we integrate by parts and use the trace estimate to see that 

'G^ < I 

r4iV-1^3 



(5.48) 



J dy-d'^G^ 




L 











< 














i?-l/2(S) 





< i|5°^u||^i/2(s) iiv: 



-47V- 1^3 



1/2 
1 1/2 



(5.49) 



i/2< ii5"«iiil|v: 

For the term we must split to two cases: ao > 1 and ao = 0. In the former case, there is at 
least one temporal derivative in so by (5.18) we have 



<{Wri\\-i/2 + <^Wn%l2)\\^ 



iN-2 



3/2 y 



dtG^ 



1 1/2 



^^2N^ ^2n'^2N + ^^2N- 



In the latter case, 9" involves only spatial derivatives; in this case we claim that 

{dy - aA^d"ri)d"G^ 



~ V^2N^2N + \/T^2N^^'iN- 



(5.50) 



(5.51) 
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Indeed, by the Leibniz rule, we have 

Q<l3<a 

For the first term, we use the integration by parts to see that 

= \ ! |5"r/|2div,u-^c7 / |9"V,r?|2div,n+ /" (75"V*r?(V*5°?7 • V*u) 

^ (ll'?ll47V-l/2 + ll'?ll4iV+3/2) ll^ll4iV+l/2 W^MIh^J^) < \/'^2N^2nK.. 

For the second term, we estimate for |/?| > 1, similarly as Lemma 5.2, 



(5.52) 



(5.53) 



1/2 



V^2V^2iV + ^-^2JV. 



(5.54) 



1/2 



Hence, we have 

In light of (5.53) and (5.55), we conclude the claim (5.51). 

Now, by (5.44)-(5.51), we deduce from (5.41) that for ah \a\ < AN with aQ<2N - 1, 

1 d 



(5.55) 



2dt 



[ |a"(xn)|'+ / \d^v\'' + <^\d''V.rj\A+l [ |D5«(xu)r 
Jn JT, / ^ Jn 



(5.56) 



The estimate (5.40) then follows from (5.56) by summing over such a. □ 
Proposition 5.7. For any e G (0, 1) there exists a constant C{e) > such that 

+ ^JV+2 ~ V ^2n'^N+2 + ^'^N+2 + C'(£)^iV+2- (5.57) 

Proof. The proof is similar to Proposition 5.4, except using (5.20) in place of (5.18) and esti- 
mating in a slight different way the term when a = 2(Ar + 2), namely, 

J {d"ri - aA^d"r])d"G'^. (5.58) 

But by (5.20), we have 

J^id^r, - aA,d^n)d'^G' < i\\d'^ri\\_y, + a \\d^ri\\,/,) || v5^+'^G'||^^^ < V^25^+2- (5-59) 

Using this, we can conclude the proposition. □ 



5.3. Comparison results. In this subsection, we shall now show that, up to some error terms, 
is comparable to and that is comparable to P„, for both n = 2N and n = N + 2. We 
begin with the result for the instantaneous energy. 



Theorem 5.8. It holds that 



and 



'^N+2 '^N+2 "T ^2N^N+2- 



(5.60) 
(5.61) 
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Proof. We first let n denote either 2N or N + 2 throughout the proof, and we compactly write 



2n = V, 



2n-2^1||2 



72rt-2^2 ||2 



2rt-2^3||2 



2n-l/^4||2 



llo+Iro lli + lr*o '-^ lli/2 + lr*o iii/2- 

(7,t 



Note that the definitions of £^ = £n + £n guarantee that 



\dMl 



3=0 



2n-2j 



j=Q 



2n-2j+l 



(5.62) 



(5.63) 



Now we let j = 0, . . . , n — 1 and then apply df to the equations in (5.14) to find 



in Q 
in Q 



{dfpl - B{dlu))e3 = {dlr, - aA^dlr,)e3 + dfG^ on S 
^^dfu = on E^. 



(5.64) 



Applying the elliptic estimates of Lemma 2.15 with r = 2n — 2_7 > 2 to the problem (5.64) and 
using (5.62)-(5.63), we obtain 



dju 



< 



2 

2n-2j 



d{p 



+ 



< 



2n-2j-2 
2 

2n-2j-3/2 
2 



2n-2j-l 
+ 



dfG' 



2n-2j-2 
2 



+ 



dlG' 



2n-2j+l/2 



+ 



2n-2j-l 



(5.65) 



2n-2j-3/2 



2n-2(i+l) 

A simple induction on (5.65) yields, by (5.63) again, 

n r, n—l r. 



j=0 



U 



2n-2j 



<\\dMl + £n+2n<£Z+2n. 



j=0 



2n-2j-l 



On the other hand, we use the boundary condition 

dtV = ^3 + on S 
to have that for j = 1, . . . , n, by the trace theory, (5.62) and (5.66), 



2n-2j+3/2 



< 



< 



^2n-2j+3/2(s) 



+ 



2n-2j+3/2 



2n-2(j-l) 



+ 



(5.66) 



(5.67) 



(5.68) 



2n-2(j-l)-l/2 



For J = n + 1, we need to use the geometric facts: d^'^^rj = D^u- N and div^(D"u) = 0. Then 
by Lemma 2.11 and (5.66), we have 



Consequently, summing (5.63), (5.66) and (5.68)-(5.69), we obtain 



CO- < OCT I 



(5.69) 



(5.70) 



Setting n = 2A^ in (5.70), and using (5.17) of Lemma 5.2 to estimate Z2N ^ (^2iv)^ + ^^2N, 
we then obtain (5.60); setting n = A?^ + 2 in (5.70), and using (5.19) of Lemma 5.2 to estimate 
2'N+2 ^ ^2Ar^w^+2' then obtain (5.61). □ 



Now we consider a similar result for the dissipation. 
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Theorem 5.9. It holds that 



and 



(5.71) 



'^N+2 ^ '^N+2 + £2n'^N+2- (5-72) 

Proof. We again let n denote either 2N or TV + 2 throughout the proof, and we compactly write 

yn 



10 



1/2 



*o 



, 2 ||Y7 2n^4||2 

1/2 + ^ Ir* <^ lll/2- 



(5.73) 



Notice that wc have not yet derived an estimate of rj in terms of the dissipation, so we can 
not apply the elliptic estimates of Lemma 2.15 as in Theorem 5.8. It is crucial to observe that 
we can get higher regularity estimates of u on the boundary S from P„. Indeed, since S is fiat, 
we may use the definition of Sobolcv norm on S, the trace theorem and Korn's inequality of 
Lemma 2.5 to obtain, by the definitions of P„ = + P+, 



din 



i?2n-2j+l/2(2) 



dtixu) 



< 



< 



j|^2n-2j + l/2(5]) 
2 



L2(S) 



+ 



//1/2(S) 

(5.74) 



This motivates us to use the elliptic estimates of Lemma 2.16. 

Now we let J = 0, . . . , n — 1 and then apply to the equations in (5.14) to find 

'-Ad{u + Vd{p = -di+\ + diG^ inn 
div dju = df in Q 

on S 
on Sfe. 



diu = o 



(5.75) 



Applying the elliptic estimates of Lemma 2.16 with r = 2n — 2j + 1 > 3 to the problem (5.75) 
and using (5.73)-(5.74), we obtain 



< 



< 



2n-2j+l 



+ 



v&Ip 



2n-2j-l 



2 

+ 


diG' 


2 

+ 


dlG' 


2 

+ 


diu 


2n-2j-l 




2n-2j-l 




2n-2i 





2n-2j+l/2 



(5.76) 



+ + 



2n-2(j+l)+l 

A simple induction on (5.76) yields, by Korn's inequality. 



j=0 



2n-2j+l 



+ 



< \\dMl+1^n+yn<'Dn + yn 



2n-2j-l 



(5.77) 



Now that we have obtained (5.77), we estimate the rest parts in We will turn to the 
boundary conditions in (5.14). First we derive estimates for rj. For the term dfrj for j > 2 we 
use the boundary condition 

dtr] = M3 + G'' on S. (5.78) 
Indeed, for j = 2, . . . , n + 1 we apply to (5.78) to see, by (5.77) and (5.73), that 



diri 



2n-2j+5/2 



< 



< 



2 

+ 

2n-2(j-l)+l 



+ 



di-^G* 



2n-2j+5/2 



(5.79) 



2n-2(j-l)+l/2 



<'Dn + yn 



(5.80) 



VISCOUS SURFACE WAVES 53 
For the term dtr), we again use (5.78), (5.77) and (5.73) to find 

\\dtr]\\ln+l/2 + \\0tV\\ln-l/2 ^ 1^3 ||^2„+i/2(s) + \\G^\\ln+l/2 + ll'^iL-1/2 

^ \H\L+l+yn<T>n+yn- 

For the remaining 77 term, i.e. those without temporal derivatives, we use the boundary condition 

-aA^r] + ri = p- 2d^u^ - G3 on E. (5.81) 

Notice that at this point we do not have any bound of p on the boundary S, but we have 
bounded in fi. Applying V* to (5.81) and then applying the standard elliptic theory, by the 
trace theory, (5.77) and (5.73), we obtain 



2 IIV7 ||2 , IIY7 ||2 

^ ll^*'7ll2n+l/2 + ll^*^ll2n-3/2 

~ II^*P|Ih2"-3/2(E) + ||V*(?3U3||^2n-3/2(-2) + 



2n-3/2 



(5.82) 



< iivpiiL-i + ikiiL+i + i|GiL-i/2 ^ + yn- 

We may then use Poincare's inequality on S (since (5.3)) to obtain from (5.82) that 

lkllL+3/2 + ll'7ll2n-l/2 ^ ll'^llo + ll^*'?ll2n+l/2 + ll^*'?llL-3/2 

^ CT^ l|V*??||2n+l/2 + l|V*?7||2n-3/2 ~ +yn. 

Summing (5.79), (5.80) and (5.83), we complete the estimates for rj: 

n+l 



(5.83) 



ll'7ll2n+3/2 + <''^ ll^t^ll2n+l/2 + ll^ll2n-l/2 + ll^t^ll2n-l/2 

j=2 



< 



2n-2j+5/2 



Vn + yn- (5.84) 



The 77 estimates (5.84) allows us to further bound d^p . Applying j = 0, . . . , n — 1 to 
(5.81) and employing the trace theory, (5.77), (5.84) and (5.73), we find 



&ip 



< 



L2(E) 



< 












2 






2 


d{r] 


+ 

2 




u 


+ 

2 



L2(E) 



<T>n + yn- 



+ 



(5.85) 



By Poincare's inequality on Q (Lemma 2.4) and (5.77) and (5.85), we have 



2 .2 .2 
dfp = dfp + V&Ip < dfp 



L2(S) 



+ 



Vdip 



<f>n + yn. 



(5.86) 



In light of (5.86), we may improve the estimate (5.77) to be 



n-l 



j=0 



U 



2n-2j+l 



2n-2j 



<'Dn + yn. 



Consequently, summing (5.84) and (5.87), we obtain 



(5.87) 



(5.88) 



Setting n = 2N in (5.88) and using (5.18) of Lemma 5.2 to estimate 3^27V ^ ^2N^2N + ^•^2iv, 
we then obtain (5.71); setting n = iV + 2 in (5.88), and using (5.20) of Lemma 5.2 to estimate 
yN+2 < ^^Ar^Ar+2' we then obtain (5.72). □ 
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5.4. Global energy estimates. In this subsection, we shall now conclude our global energy 
estimates of the solution to (1.11). We begin with the estimate of T2N- 

Proposition 5.10. There exists a universal constant Q < 5 < \ so that i/^2Ar(^) ^ then 

sup J"2iv(r) < J'2iv(0) + t [ V^j^ for allO<t<T. (5.89) 

0<r<t Jo 

Proof. Based on the transport estimate of Lemma 2.19 on the kinematic boundary condition, 
we may show as in Lemma 7.1 of [13] that 



(5.90) 



(5.92) 



sup J^2N{r) <exp 

<<r<t \ Jo 

■2iv(0) + t j\l + S^f,{r))V^r^{r)dr + (^^* ^/lC{rW^)d' 
The Sobolev and trace embeddings allow us to estimate /C < £"^^2' hence 

^/W)dr < 1^ y[£h^)dr <^6j^ __L_dr- < ^5. (5.91) 
Then by (5.91), we deduce from (5.90) that 

sup J2iv(r) < .F27v(0) +i / + sup J2iv(r) { I ^flC{r)dr\ 

0<r<t Jo 0<r<t \Jo / 

<^2Jv(0)+t / + 6 sup J-2^(r). 

JO 0<r<t 

By taking 6 small enough, we may absorb the right-hand J^2N term of (5.92) onto the left and 
deduce (5.89). □ 

Now we show the boundedness of + /q 2^2Ar- 
Proposition 5.11. There exists S > so that «/^2Ar(^) — ^' ^^^'^ 

Smit) + [ V^j^< ^2^(0) + -^2iv(0) for allO<t<T. (5.93) 
Jo 

Proof. Note first that since <£'^jv(*) — ^2Ar(^) ^ by taking S small, we obtain from (5.60) of 
Theorem 5.8 and (5.71) of Theorem 5.9 that 

£2N ^ ^ + ^•^2iv, and t>2N < ^2A' < ^2JV + ICJ'2N. (5.94) 

Now we multiply the integration in time version of (5.25) by a constant 1 + ^9 > 2 and add 
then this to the integration in time version of (5.40) to find, by (5.26), 

(1 + p)s^'^it) + s^'+{t) + [\i + ppl^ + V+j, 

Jo 

< (1 + /3)(^2t (0) + (^2^(0))=^/') + ^"2^(0) + (1 + mmit))'/' (5.95) 

+ / {2 + P) v/^25Jiv + V^fiV^-^2iV + £PfiV + <^(^)^2 
Jo 

Then by (5.94) we may improve (5.95) to be 



^2iv(i) + / 1^2N + /3^2iV < (1 + mNiO) + (1 + /3)(f2V(i))'/' + mJ'2N{t) 

+ f{2 + /?) V^^^2iV + + K.F2N + £l?2iV + C(e)^2JV (5-96) 

Jo 
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Taking e sufficiently small first, /3 sufficiently large second, and b sufficiently small third, we 
may deduce from (5.96) that 

^2N (t) + /* 2^2iv ^ £2N (0) + JCT2N + /* ^/V^^^K^ + KT^N ■ (5.97) 
Jo JO 

Since ^27v(^) — ^i is easy to use Proposition 5.10 to verify that 

/C(t) J-2Ar(i) < 5(1 + 0"^'^+^^2iV < 5J^2iv(0) + b /* V^r,{r)dr, (5.98) 

Jo 

/* }C{r)J'2N{r)dr < (5J2iv(0) + 6 f V^^{r)dr, (5.99) 
Jo Jo 
and ^ ^ 

y" ^Pf^(r)/C(r)J^2Jv(r) < J^2iv(0) + ^/ T^2N{r)dr. (5.100) 
We then plug (5.98)-(5.100) into (5.97) to have 

'?2V (i) + /* ^2^ < £2N (0) + J2iV (0) + /* P2"iv (r)c?r. (5.101) 
Jo Jo 

Taking further 6 sufficiently smaller, we may then conclude (5.93). □ 

It remains to show the decay estimates of 
Proposition 5.12. There exists 5 > so that if Q2n^) — ^> ^^^^ 

(1 + t)^^-^fXr+2(*) ^ ^2^(0) + ^2Jv(0) for all(}<t<T. (5.102) 
Proof. Since £^2Ar(*) — ^2Ar(^) — ^y taking S small, we obtain from (5.61) and (5.72) that 

^N+2 ^ ^N+2 ^ ^N+2^ P^+2 ~ ^Ar+2 ^ ^Ar+2- (5.103) 

By these estimates and the smallness of (5, we may deduce from Proposition 5.5 and Proposition 

5.7 that there exists an instantaneous energy, which is equivalent to £^^2^ simplicity is 
still denoted by £^_^_2j such that 

j^£n+2+T^n+2<0. (5.104) 
In order to get decay from (5.104), we shall now estimate <f^_|_2 in terms of T>'^_^2- Notice 

1 1 1 1 2 1 1 1 1 2 

that ^?^_|_2 can control every term in £^_^_2 except a ||??||2(iv+2)+i ll'?ll2(Ar+2)- '^^^ point 
is to use the Sobolev interpolation. Indeed, we first have that 

IkllW) ^ ll^ll2(iV+2)-i/2 hWlV^ < m+2)\£2Nf-', where 9 = (5.105) 
Similarly, we have 



ll'7ll2(iV+2)+l — ll'7ll2(W+2)+3/2 ll^ll2(Ar+2)+l/2 



(5.106) 



^ ll^ll2(Ar+2)+3/2 ll^ll2(Af+2)-l/2 ll^ll4iV ~ l^Af+2j l^27Vj 

Hence, in light of (5.105)-(5.106) we may deduce 

£%+2 < m+2)\£2Nf-'- (5.107) 

Now since by Proposition 5.11, 

sup ^2V(r) < £2n{^) + ^2iv(0) := A^o, (5.108) 

0<r-<t 

we obtain from (5.107) that 

^^^+2<-Mo"'(^^+2)'- (5.109) 
Hence by (5.104) and (5.109), there exists some constant C > such that 

|«« + ;4<^S«)'-<0. where » = i-l = j^. (5.110) 
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Solving this differential inequality directly, we obtain 

Using that £j^^2i^) ~ -^o and the fact 1/s = 4n — 8 > 1, we obtain from (5.111) that 

^N+2{t) ^ ^ ^^^^1/, ^ ^ ^1/,^ - (-^ ^ ^47V-8) ■ ^^■''^) 

This directly implies (5.102). □ 

Now we can arrive at our ultimate energy estimates for ^|jv- 
Theorem 5.13. There exists a universal < S < 1 so that i/^2iv(^) ^ ^> ^^^^^ 

G2N{t) < G2n{0) for allO<t<T. (5.113) 
Proof. The conclusion follows directly from the definition of and Propositions 5.10-5.12. □ 

5.5. The limit as a — >■ 0: the proof of Theorem 1.2. We can now prove a cr-independent 
global well-posedness of (1.11) which allow us to pass to the limit as cr ^ on the whole time 

interval [0, 00). 

Proof of Theorem 1.2. Assume that the initial data Uq and ryg satisfy the assumptions of The- 
orem 1.2 for < (7 < 1. For each cr, we let {u'^ ,p'^ ,ri'^) be the local solution to (1.11) on [0,To] 
produced by Theorem 1.1. By the estimate (5.113) of Proposition 5.13 and (4.101), we have 

Omit) < ^2^(0) < P{£S) ^ ^0 for all < t < To. (5.114) 

This allows us to employ a continuity argument as in Section 9 of [13] to prove that there exists 
a universal 5o > such that if £q < 6q, then the solution {u'^ ,p'^, rj"^) is a unique global solution 
to (1.11) on [0,00) and satisfies the a-independent global estimates (1.21). 

As for Theorem 1.1, the cr-independent estimates (1.21) yield the strong convergence of the 
whole family of [u" ,p'^ ^rj") to the unique limit (u,p,r]) which is the unique strong solution of 
(1.11) for (7 = on [0, 00) that satisfies the estimates (1.22). The strong convergence is in any 
functional space that ^2iv "^^n compactly embed into, which is more than sufficient for us to 
pass to the limit as cr —> in (1.11) for each t G [0, 00). We thus conclude the global-in-time 
zero surface tension limit and complete the proof of Theorem 1.2. □ 

5.6. Byproduct: a global well-posedness of (1 11) with surface tension. For each fixed 
0" > 0, we can establish from the proof of Theorem 1.2 a global well-posedness of (1.11) with 
exponential decay rate. In this subsection, we will allow the constants C, 60, etc. depend on a. 

Theorem 5.14. Let a > 0. Assume that the initial data uq and r]o satisfy the {2NY^ compat- 
ibility conditions (3.83), and that rjo satisfies the zero average condition (1.19). There exists a 
Sq > such that if ||ito||4jv + ll'7o|l4Ar+i < <^0; then there exists a unique global solution {u,p,rj) 
to (1.11) on [0,00). The solution {u,p,r]) obeys the estimate 

sup£2N{t)+ [ 1^2N{s)ds <C (\\uo\\Ij^ + \\r]o\\l,^,^) . (5.115) 
t>o Jo ^ ' 

Moreover, there exist constants Cx^Ci > such that 

elN{t) < CA{0)e-^^K (5.116) 
Proof. Since J^2N < C£\^, we may derive from (5.94) that 

£In < £In < S2N, and V2N < -Dl^ < 'D2N. (5.117) 
We may then deduce from Proposition 5.4 and Proposition 5.6 that there exists an instantaneous 
energy, which is equivalent to £2N^ but for simplicity is still denoted by £2^^ such that 

j/2N + 1^1n<0. (5.118) 
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Integrating the inequality (5.118) directly in time, we get (5.115). But we see 

so the decay estimate (5.115) follows by (5.118) and an application of Gronwall's inequality. □ 
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